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During  recent  years ,  logarithmically  periodic  antennas  have  been 
widely  used  due  to  their  frequency  response  characteristics,  simplicity 
of  design  and  directivity.   However,  their  theory  of  operation  still  is 
in  a  development  phase,  and  very  few  models  have  been  fully  analyzed. 
The  present  paper  is  an  attempt  to  analyze  the  operation  of  a  zig-zag 
model  that  has  the  property  of  being  symmetrical,  and  suitable  for 
operation  against  ground.   The  radiation  pattern  of  the  antenna  is 
obtained  for  different  models  of  current  distribution,  and,  finally,  the 
impedance  characteristics  and  an  approximate  current  distribution  are 
obtained,  using  non-uniform  transmission  line  theory.   The  results 
obtained  show  reasonable  agreement  with  experimental  data,  and  confirm 
conclusions  drawn  from  physical  considerations. 
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I  -  LOG-PERIODIC  ANTENNAS  -  HISTORICAL  BACKGROUND 

1.1  -  Frequency  Independent  Antennas  and  Log-Periodic  Antennas 

Frequency  independent  antennas  were  originated  in  1954,  when  V.  H. 
Rumsey  [14]  started  the  idea  that  a  structure  characterized  by  angles,  and 
without  any  characteristic  length  dimension,  should  behave  independently 
of  the  frequency  of  operation.   As  a  practical  example,  Rumsey  proposed  an 
equiangular  spiral  structure.   J.  D.  Dyson  carried  out  a  study  of  such  a 
structure  composed  of  two  plane  sheets  of  metal  extending  from  a  central 
point  according  to  this  geometry,  and  he  proved  that  the  model  not  only 
had  frequency  independence  characteristics,  but  also  that  the  interesting 
part  of  the  antenna  was  confined  to  the  area  inside  an  "active  region," 
where  the  dimensions  become  of  the  order  of  half  a  wavelength.   That  is, 
the  outward  area  is  not  necessary  and  the  structure  can  be  truncated  at  a 
convenient  size,  and  will  be  useful  within  a  range  of  frequencies  whose 
lower  limit  is  determined  by  the  peripheral  dimensions,  and  whose  upper 
limit  is  determined  by  the  spacing  at  the  feeding  terminals,  near  the 
center  of  the  antenna. 

A  structure  like  this  is  bidirectional  and  has  circular  polarization. 
Also,  since  the  active  region  rotates  along  the  spiral,  the  polarization 
rotates  too.   Strictly  speaking,  then,  not  all  the  characteristics  of  the 
antenna  are  frequency  independent,  but,  instead,  they  repeat  periodically 
every  time  the  frequency  is  changed  so  as  to  reproduce  a  reference  space 
condition.   The  input  characteristics,  pattern  shape  and  field  intensity 
are  constant. 

One  year  later,  in  1955,  R.  H.  DuHamel,  who  was  working  with  Rumsey 
in  broadband  antenna  development,  proposed  a  structure  that,  instead  of 
being  entirely  defined  by  angles,  was  interrupted  by  discontinuities 


Fig.  1  -  DuHamel's  log-periodic  teeth  structure 


properly  located.   The  first  model  built  (Fig.  1)  was  formed  by  a 
circular  sheet  of  metal,  cut  along  two  symmetrical  circular  sectors,  and 
with  teeth  placed  along  the  radial  edges.   The  size  of  the  teeth  and 
their  distances  to  the  center,  are  related  by  a  constant  ratio  p.   It  is 
evident  that  such  a  structure,  within  its  size  limits,  will  repeat  its 
radiating  characteristics  every  time  the  frequency  changes  by  p;  that  is, 
for  all  the  set  [p  f_] .   A  plot  of  these  frequencies  in  a  logarithmic 
scale  will  show  them  equally  spaced  by  log  p,  and  that  is  the  origin  of 
the  name  log-periodic  given  to  this  kind  of  geometry. 

The  model  shown  in  Fig.  1   is  a  bidirectional  antenna  with  large 
half  power  beam  width.   Bending  one  of  the  halves  of  the  antenna  toward 
the  other,  so  as  to  form  an  acute  angle,  D.  E.  Isbell  obtained  a  unidirec- 
tional antenna  with  the  surprising  characteristic  of  radiating  toward  the 
vertex;  that  is,  backward  in  relation  to  the  feed.   Meanwhile,  proceeding 
with  his  work  in  the  equiangular  spiral  antenna,  Dyson  devised  another 
unidirectional  antenna  by  wrapping  the  two  metal  sheets  of  the  planar 
spiral  around  a  cone.   Again,  he  obtained  a  backward  radiating  structure, 
whose  characteristics  could  be  adjusted  by  means  of  the  cone  angle  and 
the  spiral  rate  of  growth. 

These  two  pioneer  models  brought  log-periodic  antennas  to  the  field 
of  practical  application;  and  they,  and  their  variations,  have  been 
extensively  used  in  the  VHF  band  and  above.   They  are  particularly  suited 
for  military  applications  where  more  often  than  not,  a  single  antenna  is 
required  to  cover  a  very  large  range  of  frequencies.   As  it  might  be 
expected,  these  antennas  have  moderate  gains,  and  for  fixed  circuit 
applications  they  do  not  match  high  gain  structures.   However,  when  a  trade 
between  even  performance  along  a  large  band,  and  optimum  performance  only 


at  one  frequency  or  narrow  band  is  allowed,  the  log-periodic  models  offer 
the  best  possibilities. 

1.2  -  Unidirectional  Planar  Structures 

Traditionally,  it  is  with  the  lower  frequencies  that  the  search  for 
efficient  radiating  systems  has  been  less  successful.   Good  electrical 
performance  has  been  obtained  at  the  expense  of  size,  and,  although 
theoretically  it  would  be  possible  to  build  optimum  systems,  the  costs 
involved  in  the  ground  occupied  and  the  hardward  employed  make  such  attempts 
prohibitive.   In  the  HF  band,  the  best  compromise  between  size  and  opera- 
tional characteristics,  during  the  past  three  decades,  has  been  the  rhombic 
antenna.   Today,  the  rhombic  still  is  the  practical  structure  with  highest 
gain,  and  when  used  with  a  range  of  frequencies  within  a  ratio  of  not  more 
than  3:1,  it  still  might  be  the  best  answer  for  an  HF  radio  link.   Its 
main  inconveniences  are  the  size  and  the  variable  take-off  angle,  but  the 
latter  can  be  used  with  advantage  by  appropriate  selection  of  frequency 
according  to  ionospheric  conditions. 

The  log-periodic  models  of  the  type  mentioned  above  are  not,  in 
general,  feasible  for  HF  long  range  communications.   Besides  being  too 
bulky,  they  employ  sheets  of  metal  that  would  make  the  cost  too  high. 
Two  important  developments  were  to  come,  and  solve  these  two  problems. 

The  first  was  when  DuHamel  and  his  research  team  verified  that  only 
the  peripheral  strips  of  the  metal  sheets  are  active  in  the  radiation 
mechanism,  so  that,  essentially,  all  the  structures  can  be  reduced  to 
wire  frames. 

The  second  was  when  Isbell  [11]  built  the  first  array  of  log-periodic 
half-wave  dipoles  (Fig.  2a).   Again,  here  a  surprising  fact  was  verified. 
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Fig.  2a  -  Isbell's  transverse  dipole  array 
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Fig.    2b  -  A  tapered-ladder   for   operation  above  ground 
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The  structure,  as  all  log-periodic  structures  so  far  realized,  radiates 
backward;  but,  to  obtain  proper  operation,  every  other  dipole  of  the 
array  is  fed  through  a  reversal  of  the  feed  line  connections. 

This  dipole  array  established  the  applicability  of  log-periodics  in 
the  HF  band.   The  structure  can  be  built  for  vertical  or  horizontal 
polarization,  and  in  the  latter  case,  by  adequate  choice  of  the  slope  of 
the  feed  line,  either  fixed  or  variable  take-off  angles  can  be  obtained. 
The  first  models  had  tubular  dipoles,  with  radius  scaled  to  the  structure 
period.   Granger  Associates,  who  might  be  properly  recognized  as  the 
pioneers  in  large  scale  building  of  HF  log-periodic  antennas,  have 
replaced  the  tubular  dipoles  by  wire  triangular  elements  with  appreciable 
reduction  in  weight,  making  it  possible  to  build  the  rotatable  long  range 
log-periodic.   This  is  a  steerable  antenna,  mounted  on  a  pole  high  enough 
so  that  the  take-off  angles  are  as  low  as  10°,  and  light  enough  to  be 
rotated  by  a  small  electric  motor  through  a  360°  azimuth  span. 

Incidentally,  one  of  the  problems  that  arose  with  these  antennas  was 
the  matching  of  the  feeding  line.   The  antenna  terminals  show  a  balanced 
input  impedance  of  the  order  of  200  ohms,  and,  typically,  a  maximum  VSWR 
with  respect  to  50  ft  below  2:1  through  the  band  of  operation.   Two  hundred 
ohms  is  not  the  most  convenient  value  to  carry  energy  from  a  distant 
transmitter,  and  the  problem  was  that  of  obtaining  a  matching  device  that 
would  enable  the  use  of  50  ohms  standard  coaxial  cable.   That  need  speeded 
up  the  research  in  the  area  of  baluns  and  matching  transformers  with 
outstanding  results.   Granger  Associates  have  supported  their  achievements 
in  antenna  production,  supplying  also  baluns  with  the  required  broadband 
characteristics,  and  multicouplers  that  enable  the  simultaneous  use  of  one 
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single  antenna  with  a  number  of  independent  transmitters;  the  total  power 
being  the  only  limitation  in  the  system  plus  a  small  separation  in  the 
frequencies  of  operation. 

Looking  for  more  and  better,  the  next  step  in  log-periodics  was  to 
try  to  obtain  a  structure  that  would  use  a  ground  image  with  the  purpose 
of  reducing  its  size.   The  reversal  in  dipole  connections  excluded  that 
possibility  in  Isbell's  array,  except  by  the  use  of  inverting  transformers 
between  radiators,  which  does  not  seem  to  be  a  very  practical  idea.   The 
more  significant  successful  structures  are  the  bent  log-periodic  zig-zag, 
the  tapered  ladder,  and  the  monopole  array.   The  bent  zig-zag,  developed 
by  J.  Greiser  and  P.  Mayes  [8],  [9],  will  be  the  main  subject  of  this 
paper.   The  tapered  ladder,  developed  by  Wickersham,  consists  of  an  array 
of  grounded  vertical  radiators  (Fig.  2b)  capacitively  coupled  to  an 
inclined  transmission  line,  the  slope  of  which  is  adjusted  for  proper  feed 
of  each  radiator.   A  variation  of  this  arrangement,  that  uses  triangular 
wire  radiators,  has  been  commercialized  by  Granger  Associates.   The 
monopole  array,  developed  by  D.  Berry  and  F.  R.  Ore  [23]  uses  phasing 
stubs  between  every  two  radiators,  scaled  with  the  structure  period.   The 
stubs  are  horizontal  and  open-circuited,  so  that  only  one  wire  is  actually 
necessary,  since  its  image  provides  the  other,  and  the  radiators  are 
vertical  tubes,  of  a  quarter  of  a  wavelength,  for  the  respective  operating 
frequency. 

In  general,  the  half-size  structures  do  not  have  as  good  a  performance 
as  the  transverse  dipole,  and  the  recent  advances  in  construction 
techniques  have  made  the  latter  a  stronger  competitor,  even  for  frequencies 
as  low  as  2  MHz, 
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1.3  -  Developments  in  the  Theory  of  Log-Periodic  Antennas 

The  development  of  log-periodic  structures  during  their  first  years 
was  mainly  the  result  of  good  intuition  and  extensive  laboratory  testing 
of  a  large  number  of  proposed  models.   As  a  rule,  only  after  a  successful 
model  was  obtained,  was  a  theoretical  foundation  to  explain  its  perform- 
ance developed,  and  only  relatively  recently  models  based  on  sound 
theoretical  investigation  are  being  realized. 

There  have  been  three  main  approaches  to  the  problem  of  analysis  of 
operation  of  LP's.   The  first,  which  we  can  call  the  classical  approach, 
uses  the  theories  developed  for  the  components  of  the  antennas  and  combines 
them  to  suit  the  particular  model  under  consideration.   The  outstanding 
example  of  this  approach  is  Carrell's  analysis  of  the  transverse  dipole 
array  [10] ,  who  described  the  behavior  of  the  structure  in  terms  of 
currents,  voltages,  input  impedance  and  far-field  with  very  close  agree- 
ment with  experiments.   This  approach  usually  has  the  limitation  of  being 
valid  for  a  particular  structure  only. 

Another  approach,  that  has  been  pursued  by  DuHamel   [15],  [16], 
Armstrong  [16],  Mittra,  Jones  [17],  [18],  [19],  and  others  is  to  develop 
an  adequate  body  of  information  and  knowledge  in  log-periodic  networks, 
on  which  a  general  theory  of  the  antenna  problem  could  be  based ■   Great 
progress  has  been  made  toward  this  purpose,  and  results  obtained  have 
helped  in  a  better  understanding  of  log-periodics,  and,  in  some  cases, 
provided  new  ideas  for  future  realizations. 

The  third  approach,  followed  by  Mayes,  Deschamps,  Patton  [20],  Greiser 
[8],  [9],  Oliner,  Mittra,  Jones  [22]  and  others  is  an  application  of 
concepts  developed  for  surface  and  leaky  wave  structures,  that  are 
typically  uniformly  periodic  structures,  to  the  active  region  of  log-periodics 
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that  can  be  approximated  as  a  uniformly  periodic  region  for  limited 
purposes.   Since  such  concepts  give  an  excellent  insight  into  the  main 
operational  characteristics  of  these  structures,  we  will  briefly  review 
some  of  its  basic  points  [4],  [5],  [7]. 

Consider  a  structure,  such  as  a  waveguide,  with  a  periodic  grating 
in  one  face,  and  suppose  that  the  waveguide  carries  an  electromagnetic 
wave.   Due  to  the  discontinuities  in  the  grated  face  there  will  be 
leakage  of  energy  through  them,  and,  also,  the  fields  inside  the  guide 
will  adjust  to  this  periodic  configuration.   According  to  Floquet's 
theorem,  "for  a  given  mode  of  propagation  at  a  given  steady-state 
frequency,  the  fields  at  one  cross  section  differ  from  those  one  period 
away  only  by  a  complex  constant"  [7].   This  can  be  expressed  as 

where  L  is  the  period  of  the  grating  and  y  is  the  propagation  factor. 

Now,  if  the  longitudinal  dimension  of  the  guide  is  made  to  coincide 
with  the  z  axis,  and  we  expand  the  equation  in  a  Fourier  series,  we  get 

E  =  E(x,y,z)e-YZ  =   £   E  (X,y)e^  (2™/L)z  e^Z,  (1) 

all  n  n 

that  is,  we  obtain  a  set  of  n  space  harmonics,  each  one  having  propagation 
constant 

.  2-rm 

r  =  y  +  j  —  . 

If  we  now  assume  that  the  system  is  lossless,  we  obtain 

where  3n  is  the  medium  phase  constant.   Now,  we  apply  this  result  to  the 
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usual  description  of  waves  bounded  by  guides.   In  general,  the  fields 
are  described  by  the  wave  numbers  and  a  relation  of  the  form 

2    2    2    2 
K  =  K  +  K  +  K   . 
x    y    z 

When  the  x  =  0  plane  is  a  ground  plane,  that  is  the  general  case 
of  interest  for  the  present  purpose,  the  y  number  is  zero,  and,  with  the 
assumption  of  no  losses,  we  have 

2    2    2    2 

B  *  K  +  3  +  B 
x    z 

where  M  K  is  the  free  space  phase  constant. 

If  now  we  recognize  that  the  x  and  z  phase  constants  are  given,  for 
the  grated  guide,  by  a  series  of  the  form  of  equation  (2),  we  obtain  the 
general  form  of  the  solution.   But,  before,  let  us  consider  an  angle  6, 
as  being  the  angle  between  the  z  axis  and  the  direction  of  a  z  component 
wavefront.   Then  we  can  relate  the  x  and  z  wave  numbers,  or  the  phase 
factors  as 


B  =  K  sin  G 
x 


and 

B  =  K  cos  6   , 

so  that  the  direction  of  leaking  energy,  that  is  the  direction  of  the 
wavefront,  is  given  by 

cos  8  =  J  •  (3) 
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Now,  according  to  (2)  we  recall  that  6   can  take  any  discrete  value 

z 

6   (note  that  n  ranges  from  all  negative  values  to  all  positive  values), 
and,  being  so,  the  ratio  (3)  can  take  an  equal  number  of  values,  positive 
or  negative.   When  this  ratio  falls  in  the  range  -1  to  +1,  that  corresponds 
to  a  solution  of  the  cos  function,  and  also  to  real  solutions  of  space 
harmonics  for  which  there  is  leaking  energy.   In  particular,  when  (3)  takes 
the  value  -1,  energy  is  fired  backward,  and  this  corresponds  to  the  mode 
of  interest  in  the  present  analysis. 

Since  the  amplitude  of  the  harmonics  decreases  with  increasing  order, 
we  would  like  to  have  the  structure  operating  in  a  low  order  mode.   For 
the  n  =  -1  mode  we  have 


6-l   ~   60   "     L 

and 

n            60  -   2tt/L 
cos  e_1  =          K 

In  log-periodic  structures,  always  there  is  some  sort  of  guide  that 
carries  energy  from  the  feed  point  toward  the  radiators,  and  discontinuities 
placed  along  the  guide  force  energy  to  leak  away.   In  the  transverse  dipole 
the  line  is  the  feeder,  a  transmission  line,  and  the  grating  is  represented 
by  the  dipoles;  in  the  teeth  structure  of  Fig.  1,  the  plane  sheets  form 
the  guide  and  the  teeth  form  the  grating;  the  bent  zig-zag  is  a  transmission 
line  deformed  by  the  radiating  elements.   In  each  of  these  models,  and,  in 
general,  in  each  log-periodic  model,  we  can  consider  three  distinct 
regions:   the  active  region  comprising  the  two  to  five  radiating  elements 
near  resonance  for  the  operating  frequency,  the  feeding  region,  between 
this  and  the  source  where  radiation  is  small  and  in  general  the  radiators 
can  be  considered  as  a  capacitive  loading,  and  the  end  region,  that  we 
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assume  is  not  reached  by  any  significant  amount  of  energy,  since  most 
of  it  was  radiated  by  the  active  region.   Since  the  active  region  is 
small,  we  now  assume  that  it  is  a  uniformly  periodic  structure,  and  then 
apply  it  to  the  above  results  for  leaky  waves. 

Perhaps  a  clearer  picture  of  such  an  application  is  by  means  of 
dispersion  diagrams.   Fig.  3  shows  a  dispersion  diagram  for  equation  (3), 
assuming  that  6n/K  >  1.   It  is  a  well-known  fact  in  transmission  theory 
that  propagation  within  boundaries  takes  place  when  the  phase  constant 
in  the  medium  is  larger  than  the  phase  constant  in  the  boundary  material, 
i.e.,  when  the  propagation  velocity  in  the  medium  is  lower  than  the 
propagation  velocity  in  the  boundary  material.   For  metallic  boundaries 
the  velocity  is  usually  taken  as  that  of  free  space,  and  we  represent 
that  condition  by  the  phase  constant  K.   This  situation  is  usually 
identified  by  the  designation  of  slow  waves  and  it  is  said  that  periodic 
structures  have  stop  bands  for  fast  waves.   These  fast  waves  that  do  not 
propagate  within  the  structure  are  the  ones  that  leak  away  and  form  the 
radiated  field. 

In  Fig.  3  we  represent  the  fundamental  and  the  first  and  second 
negative  harmonics  as  given  by  equation  (2) ,  and  we  show  the  +45°  slope 
lines  that  separate  zones  of  slow  waves  from  zones  of  fast  waves.   We 
see  that  the  fundamental  is  a  slow  wave  for  all  frequencies,  and  so  are 
all  the  positive  order  harmonics.   The  negative  order  harmonics  are 
slow  waves  for  low  frequencies,  but  as  the  frequency  increases  to  near 
the  point  B,  the  first  negative  harmonic  approaches  the  negative  slope 
line.   On  the  line  we  have 

6_1  60-2u/L 

cos  e_^  =  — =-■ — - —  =  -l  , 
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Fig.  3  -  Dispersion  diagram 
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that  is,  we  start  having  radiation  backwards;  and  until  the  positive 
slope  line  is  reached  the  radiation  angle  rotates  from  180°  to  0°. 
Meanwhile,  the  other  negative  order  harmonics  were  having  similar 
behavior  in  the  region  that  we  can  visualize  to  the  left  of  the  diagram. 

Since  the  analysis  is  applied  to  a  uniformly  periodic  structure, 
we  only  have  well-defined  active  bands,  the  first  of  which  is  the 
strongest,  separated  by  stop  bands  (these  have  the  opposite  meaning  of 
the  stop  bands  of  transmission  theory).   However,  in  log-periodic 
antennas  the  structure  is  scaled  almost  continuously  (typical  scale 
factors  range  from  .80  to  .98),  so  that  there  is  a  continuous  active 
band,  limited  only  by  the  operating  frequencies  of  the  first  and  the  last 
of  its  elements.   The  frequency  independence  characteristics  of  the 
antenna  depend  on  how  well  its  discrete  scaling  approaches  a  continuous 
scaling,  and  the  over-all  antenna  performance  can  be  predicted  from  this 
type  of  analysis.   Greiser  and  Mayes  [8],  [9]  have  used  this  technique 
with  remarkable  success,  in  predicting  whether  or  not  different  models  of 
bent  zig-zag  antennas  would  radiate,  and  how  the  radiation  characteristics 
could  be  improved.   Mittra  [17]  has  conducted  a  similar  analysis  for  the 
transverse  dipole,  showing  by  this  method  how  the  phase  reversal  was 
necessary  to  make  the  dispersion  diagrams  approach  the  radiation  zone. 

Again,  there  is  a  remarkable  lack  of  dispersion  diagrams  developed 
for  the  type  of  structures  used  in  log-periodics,  and,  as  yet,  the  theory 
was  not  really  adapted  for  them,  except  under  the  heuristic  approach 
described  above.   We  should  note  here  that,  although  this  application  of 
leaky  waves  theory  in  a  wired  log-periodic  structure  might,  at  first 
sight,  seem  rather  displaced,  there  is  reason  to  believe  that  the  theory 
is  actually  valid.   R.  L.  Bell,  C.  T.  Elfving  and  R.  E.  Franks  [21] 
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reported  measurements  of  fields  around  a  log-periodic  structure  that 
show  the  existence  of  a  forward  wave  progressing  from  the  source  along 
the  structure,  and  a  backward  wave  with  a  90  degrees  spacial  phase 
difference.   The  forward  wave  is  supported  by  the  structure;  the  backward 
wave  is  launched  by  it  into  space.   That  is  exactly  the  type  of  leaky 
wave  behavior  used  in  the  above  theory. 

The  mathematical  tools  developed  to  deal  with  the  various  problems 
posed  by  electromagnetic  theory,  although  having  in  common  the  basic 
source  given  by  Maxwell's  equations,  follow  different  courses  and  tend 
to  place  in  distinct  categories,  problems  that  have  much  in  common.   It 
seems  to  be  a  marked  trend  in  recent  years,  to  reunite  electromagnetic 
theory  in  the  sense  of  extending  concepts  and  ideas  developed  for 
particular  systems  to  other,  apparently  different,  applications. 

As  a  final  remark,  we  note  that,  since  log-periodics  are  scaled  up 
from  the  feeder,  the  active  elements  are  the  ones  in  the  active  region 
together  with  the  smaller  ones.   Smaller  elements  in  all  cases  tend  to 
increase  the  propagation  velocity  (the  feeder  approaches  a  uniform  line 
toward  the  feed  point).   If  we  would  attempt  to  build  a  log-periodic  for 
forward  radiation  we  would  scale  the  structure  in  the  opposite  sense, 
and  would  count  on  positive  order  waves  for  the  radiation  process. 
However,  while  in  the  usual  scaling  the  positive  harmonics  find  a  high 
impedance  in  the  first  part  of  the  feeder  (they  are  in  the  stop  band  of 
transmission  theory  and  mismatched  to  the  available  radiators),  in  the 
opposite  scaling  the  negative  harmonics  would  still  reach  the  end  part  of 
the  feeder  and  cause  a  disturbing  radiation  from  that  region. 
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II  -  THE  BENT  LOG-PERIODIC  ZIG-ZAG  ANTENNA  (BLPZZ) 

The  BLPZZ  (Fig.  A-l)  was  developed  by  Greiser  and  Mayes  using  the 
concepts  of  dispersion  diagrams.   As  mentioned  before,  this  was  one  of 
the  first  attempts  to  obtain  a  LP  antenna  with  symmetrical  characteristics 
that  could  be  fed  against  ground.   For  that  reason,  it  seems  particularly 
suited  for  HF  applications. 

The  basic  geometry  of  the  antenna  is  derived  from  a  zig-zag  antenna 
bent  along  a  convenient  line,  so  as  to  obtain  adequate  radiating  elements 
and  proper  connecting  stubs.   Since  the  structure  is  progressively 
elevated  from  ground  the  connecting  stubs  are  actually  radiating,  but, 
for  small  elevation  angles  that  can  be  neglected.   Their  essential  function 
is  to  provide  the  proper  phase  delay  between  radiators,  and  they  can  be, 
and  have  been,  replaced  by  other  forms  of  delay  line,  or  even  removed, 
without  destroying  the  basic  characteristics  of  operation  of  the  antenna. 

The  antenna  is  usually  defined  by  the  radiators  spreading  angle  ol, 

Ci 

the  stubs  spreading  angle  a  ,  the  elevation  angle  e,  the  period  p,  and 

s 

the  length  of  a  wire,  2b,  or  the  distance  of  one  element  to  the  geometric 
origin.   In  Appendix  A  we  develop  the  relations  among  these  parameters  and 
other  useful  quantities.   We  note  that  the  spreading  angles  and  the  period 
are  interrelated  through  the  slope  of  the  wires.   Although  each  one  of 
them  can  be  independently  changed,  that  involves  a  change  in  the  dimensions 
and  inclination  of  the  wires.   Also,  we  note  that  each  element,  for  small 
elevation  angles,  can  be  compared  to  a  rhombic  antenna  made  of  small 
wires,  and  with  a  very  large  apex  angle. 

This  antenna  has  the  interesting  characteristic  of  being  fed  in 
series,  i.e.,  the  radiators  are  part  of  the  feed  line.   The  other  planar 
LP  structure  so  far  analyzed — the  transverse  dipole — has  the  radiators 
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fed  in  parallel,  and  it  might  be  worthwhile  to  compare  the  results 
obtained.   On  the  other  side,  this  one,  also  is  a  structure  that  seems 
adequate  to  analyze  with  conventional  tools. 

In  the  first  part  of  this  work  we  obtain  a  set  of  equations  for  the 
radiation  pattern  of  the  antenna.   Then  we  assign  to  the  wires  different 
types  of  current  distribution,  and  compare  the  results  with  the 
experimental  data  obtained  by  Greiser  and  Mayes. 

In  the  last  part  of  the  work  we  consider  the  antenna  and  its  image 
as  a  non-uniform  transmission  line  and  we  obtain  an  approximation  for 
the  input  impedance  and  current  distribution. 

2.1  -  The  Radiation  Pattern  Equations  of  the  BLPZZ 

To  obtain  the  radiation  pattern  of  the  BLPZZ  structure  above  ground 
we  begin  by  considering  one  element  at  a  time.   We  define  as  an  element, 
one  period  of  the  structure.   Since  the  elements  are  related  by  a  constant 
size  factor,  that  is,  the  period  p,  the  general  equations  obtained  for 
one  element  are  valid  for  the  other  elements  as  well. 

In  each  element,  we  assume  that  the  connecting  stub  does  not  radiate, 
and  we  split  the  radiating  section  in  two  straight  wires.   Then  we 
obtain  the  fields  due  to  these  two  wires  and  their  images,  all  referred 
to  an  element  phase  center,  and  add  them  vector ially^ 

We  make  two  more  assumptions.   First,  we  assume  that  there  is  no 
mutual  coupling  between  elements;  second,  we  assume  that  the  wires  carry 
travelling  waves.   This  corresponds  to  considering  the  structure  and  its 
image  as  forming  a  radiating  transmission  line  properly  terminated;  that 
is,  without  standing  waves,   Later  we  will  further  elaborate  on  this 
point. 
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After  obtaining  the  element's  equations  we  establish  an  array 
expression,  accounting  for  spacings  and  feed  currents. 

We  outline  here  the  development  of  the  equations  for  one  wire,  wire 
1  in  Fig,  4.   Details,  as  well  as  the  work  concerning  other  wires,  can 
be  found  in  Appendix  A, 

In  general,  the  distant  field  due  to  the  flow  of  current  along  a 
wire  is  given  by 

_      .    0  -jwt  r     sin  6   -jkr  ,   v  ,  , , , 

E  ,  =  -J  T~  e     J   — e     y(s)  ds  (4) 

6  0 

where  9',  s,  1  and  r  are  defined  in  Fig.  4,  and  w,  k,  y  and  tt  have  the 
usual  significance.   y(s)  is  the  current  distribution.   In  general,  this 
will  be  given  by 

y(s)  =  IQ  e  -(a+^)s  (5) 

where  a  is  the  attenuation  factor  and  $  is  the  propagation  factor,  along 
the  wire. 

For  large  distances  the  angle  6  and  the  distance  r  are  essentially 
independent  of  the  position  along  the  wire,  so  that  we  can  consider  them 
as  constants  in  regard  to  s,  provided  that  we  account  for  the  phase  delay 
at  the  far  field  due  to  the  different  points  along  the  wire, 

Again,  referring  to  Fig.  4,  if  we  consider  the  distance  R  from  the 
center  of  the  wire,  we  can  express  r  as 

r  =  R  -  s  cos  6' 
and,  replacing  in  equation  (4),  and  defining 

I0  -  M(-  f> 

we  obtain 
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Fig.  5  -  Geometrical  relations  in  the  BLPZZ 
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Now  we  apply  this  equation  to  wire  1  in  Fig.  5.   Here  1  =  2b,  and,  then 


f  e  ds    -    e -  e  

J  -  [K  +  dk(|-c0sfl;j] 


Defining 


wyn        -j(kR+wt)  -j(kR+wt) 

F  ■  "  JT  f  "  -  60  f 


we  obtain 


_e 


Mi£I.e         ^__  ...     -«.,„  (7) 


Since  we  are  going  to  use  this  equation  with  a  digital  computer,  this 
is  a  convenient  form. 

Now  we  need  to  convert  this  equation  to  a  system  of  coordinates  to 
be  used  with  all  the  structure.   We  choose  spherical  coordinates,  as 
shown  in  Fig.  6,  and  we  obtain  the  following  transformation  equations  (see 
Appendix  A,  Parts  II  and  III): 
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cos  6'  =  cos  (y-e)  cos  0  -  sin  (y-e)  sin  9  cos  0  (8) 


,  cos  (y-e)  sin  8  +  sin  (y-e)  cos  8  cos  0  ,    . 

)1  "  ^1  sin9^  cy) 


sin  (y-e)  sin  0 
h*01   fc9j_     sin  e'±  UU; 

At  this  point  it  is  convenient  to  introduce  the  correction  to  the 
"phase  centers"  of  the  element.   Here  the  designation  "phase  center"  is 
used  to  mean  the  reference  point  of  the  element.   The  actual  phase  center 
was  not  computed,  although  physical  considerations  would  lead  to  the 
point  C  in  Fig.  5,  that  was  chosen. 

Again,  we  assume  that  the  angles  0  and  0,  as  well  as  the  distance  R 
are  not  altered  to  an  appreciable  degree,  and  we  only  take  into  account 
the  phase  delay  in  the  far  field.   The  distance  c.  ,  from  the  center  of 
the  wire  to  point  C  is  given  by  (see  Appendix  A,  Part  I) 

m   b  sin  (ye) 
1     cos  6, 


1 


where 


.       -1  L  sin  e  +  b  cos  (y-e) 

6.  =  tan ; -. r-*-1 u 

1  b  sin  (y-e) 

The  phase  difference  in  the  far  field  is  given  by  (see  Fig.  5) 


R1  -  R  =  c.  cos  a, 


where 


cos  a,  =  -  [sin  61  cos  8  +  cos  6.  sin  8  cos  0]  , 
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Fig.  6  -  Coordinate  system 


Fig.  7  -  Correction  for  reference  point 
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so 


^u   ^   ^u      c  ~jkR  .      ..    ,,v  .       ,        -ik(R+cn cosa, ) 

that  the  factor  e     in  equation  (6)  becomes  becomes  e       1    1  , 

and,  finally,  combining  equations  (7),  (8),  (9),  (10),  and  this  correc- 
tion, we  obtain 


^(l.cco,*,)  HK£-"S<#  ■W'G-10*)) 

£A  =  i  k£  L  e  _e =_e 


L 


V<*  ; 


rf[rtstt-0  $«»♦»»  (K-0"sflcosd       ^' 


and 


Now  combining  these  equations  with  the  similar  equations  obtained 
for  the  other  three  wires  (see  Appendix  B)  we  have,  for  each  element, 


n     1     2     3     4 


\  "  \  +  \  +  \  +  \  <Wb) 


For  the  whole  structure  we  will  have 

E6=  I      anE9  (15a) 

n=l      n 


N 

Ert  =    I  a  E«l                                                (15b) 

0    Ln  n  0 

n=l  n 
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where  a  is  a  factor  accounting  for  the  spacing  and  feed  currents,  that 
can  be  split  in  the  product  of  two  terms,  one  accounting  for  each. 

The  feed  current  for  the  n+1  element  is  the  feed  current  for  the 
order  n  element  delay  attenuated  by  the  length  of  wires  1  and  2  and  the 
connecting  stub.   As  mentioned  before,  the  connecting  stub  is  assumed 
to  be  a  uniform  transmission  line,  so  that  it  is  only  responsible  for  an 
additional  delay. 

Referring  to  Fig.  5,  the  distance  along  the  wires  from  the  feed 
point  of  element  1  to  the  feed  point  of  element  2  is 


2b1  +  2b1  +  2b^  +  2b  ^  =  4bx  +  4b  ^  =  4(b1  +  bp 


From  Appendix  A  we  have 


,  ,   _1   tan  as 
1    y—  tan  aE 


so  that  the  distance  becomes 


d2  ■  4bi  r + 


tan  as 
/p  tan  aE 


The  distance  to  the  feed  point  of  element  3  is 

+  4b„ 


d3  =  4b1 


tan  as 


/p  tan  aE  i 


1  + 


tan  as 

/p  tan  aE  - 


4b. 


1  + 


tan  as 
/p  tan  aE 


i  +  i 

P 


,  since  b„  =  b. ; 


and  the  distance  to  the  order  n  element  is 


d  =  4b, 
n     1 


1  + 


tan  as 
/p  tan  aE 


1  +  -  +  ^7  +  . 
P    2 
P 


.  + 


n-2 


=  4b    1  +  tan  as 


SUM 


/p  tan  aE 
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where  SUM  =  1  +  -  +  -^  +  .  .  .  + 


p    2  '  ••*  '   n-2  * 
P         P 


If  we  subtrace  —  from  S,  where, 
P 


p   p    2    3         n-1 
P    P         P 


we  get 

S(i  -I)  =  i  -I—   As  =  1-  ^f"1 
p         n-1  1  -  1/p 


and   the  distance  becomes 

d     -  4b,    [  1  ■+  -SSSLiSS I    1"  1{Pn"1 

1  I     /7  tan  aE  I   1  "  ^ 

If  we  call  3  the  propagation  constant  along  the  radiating  wires,  and 
assume  that  along  the  stubs  the  constant  is  K,  as  in  free  space,  we  have 
the  phase  angle  of  the  feeding  current  given  by 

4biK||  +  ^2^)ifL^i  (16) 

»     vp  tan  aE  ' 

We  could  obtain  an  analogous  expression  for  the  attenuation,  assuming 
that  the  attenuation  factor  was  constant  along  the  radiating  structure; 
however,  physical  considerations  imply  different  attenuations  along  the 
line,  so  that  we  obtain,  for  the  order  n  element  feed  current  an  attenua- 
tion given  by 

n-1 

y   4b. a.  (17) 

i-1    X  X 

where  a.  is  the  average  attenuation  for  each  element. 

Then  the  correction  factor  for  current  in  the  array  expression  is 
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Fig.  8  -  Correction  for  array  spacing 
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To  introduce  now  a  correction  for  spacing,  we  are  going  to  use  as 
the  reference  point  the  "phase  center"  of  the  array.   Again,  we  do  not 
know  exactly  where  the  phase  center  is,  but  we  expect  in  that  respect  a 
behavior  similar  to  other  log  periodic  structures  [10],  [21].   Being  so, 
the  phase  center  should  lie  by  the  element  that  is  nearest  to  resonance, 
and  that  we  call  the  M  element.   Then,  according  to  Fig.  8,  we  obtain 


R  =  R  -  (L  cos  e  -  L  cose)  cos  9" 
n         M         n 


where  cos  6"  =  sin  6  cos  0  and 

R  =  R  -  (L.„  -  L  )  sin  8  cos  0  cos  e  (19) 

n         M    n 

...        -ikR  ,        -ikR  +ik(L.  -L  )sin  0  cos  0  cos  e 
In  equation  (6) ,  e  J   becomes  e     e     w  n 

and  we  insert  in  the  array  expression  the  factor 

jk(L-L  )sin  9  cos  0  cos  e  or\\ 

Then  the  total  correction  a  is 

n 


Finally,  the  array  equations  become 
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(2  2a) 


*-£■  «E>» 


(2  2  b) 


And  from  here  we  obtain  the  total  field  E  by  taking  the  real  parts  of  E 
and  E^,  and  combining  them  vectorially,  so  that 


=  IEIZ1 


where  V   is  an  angle  counted  positive  in  the  clockwise  sense  starting  in 
the  direction  of  the  a_,  vector  by  an  observer  at  the  array  phase  center. 

2.2  -  Current  Distribution  -  Radiation  Diagrams 

Before  programming  the  equations  obtained  in  section  2.1  in  a 
computer,  we  had  to  decide  on  the  current  distribution  that  would  be 
most  closely  approximate  to  the  real  situation.   With  that  purpose  in 
mind,  we  made  the  following  considerations: 
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Elememt" 


Note:   n  denotes  resonant  element 


Fig.  9  -  Models  of  current  distribution 

(1)  As  obtained  in  other  LP  networks 

(2)  Assuming  travelling  wave  attenuated  at 
active  region 
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(1)  The  antenna  and  its  image  form  a  radiating  transmission 
line,  and  since  experiments  [8],  [9]  show  that  its  input  impedance 
is  mainly  resistive  and  almost  constant  along  the  operating  band, 
having  a  small  capacitive  reactance  probably  due  to  input  termination 
effect,  it  is  reasonable  to  assume  that  the  line  is  well  terminated. 

(2)  The  elements  forming  the  feed  region  must  see  the  same 
impedance  toward  the  end,  and  that  is  the  input  impedance  of  the 
antenna.   This  is  so  because  that  impedance  does  not  essentially 
change  with  frequency,  and  a  change  in  frequency  is  equivalent  to 

a  change  of  antenna.   That  is,  when  we  decrease  the  frequency  by  one 
period,  we  move  the  active  region  one  period  toward  the  end,  and 
that  is  equivalent  to  using  one  antenna  "one  period  larger;"  the 
impedance  that  was  seen  by  the  input  terminals  must  now  be  seen  at 
the  end  of  the  first  element.   But  since  the  input  impedance  remains 
the  same,  that  means  the  first  element  is  a  good  match,  and  so  must 
be  all  the  elements  forming  the  feed  region.   Then  the  current  along 
the  feed  region  must  be  essentially  constant,  except  for  a  relatively 
small  amount  of  attenuation. 

(3)  The  energy  reaching  the  end  region  must  be  very  small.   The 
structure  is  abruptly  terminated  without  any  dissipation  load,  and 
it  has  been  shown  experimentally  [8] ,  [9]  that  the  effect  on  the 
radiation  pattern  of  ending  the  structure  by  a  short  or  an  open 
circuit  is  very  small,  for  frequencies  well  above  the  lowest  frequency 
of  operation.   Then,  it  seems  reasonable  to  assume  that  the  current 
dies  off  after  the  active  zone,  and  is  essentially  zero  at  the  last 
elements. 
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(4)  The  current  distribution  at  the  active  zone  is  more 
problematic.   The  non-uniform  line  forming  the  antenna  is  continuous 
in  the  sense  that  its  impedance  [2]  changes  without  interruption, 
and  has  the  property  of  a  constant  product  of  reactance  times 
susceptance.   We  can  call  this  a  continuously  tapered  line,  with 
linear  taper,  of  the  type  often  used  to  match  impedances.   Under 
this  point  of  view  it  is  continuously  matching  itself. 

On  the  other  side,  if  we  think  of  the  line  as  composed  of  series 
sections  connected  by  stubs,  then  we  have  to  consider  the  fact  that 
some  of  these  sections  (the  ones  forming  the  active  region)  are  near 
their  natural  resonance  frequency  and  are  likely  to  form  standing  waves, 

The  first  point  of  view  would  lead  us  to  a  heavily  attenuated 
travelling  wave.   The  second  leads  us  to  the  type  of  current  distribu- 
tion found  by  Carrel  [10]  in  the  dipole  array,  and  by  Mittra  and 
Jones  [18]  for  a  theoretical  line.   A  point  to  note  here  is  that  both 
these  cases  had  parallel  loads,  while  the  present  case  is  of  series 
loading.   Fig.  9  shows  the  two  types  of  current  distribution  corre- 
sponding to  the  two  points  of  view. 

We  started  by  trying  these  two  models.   The  results  obtained  can  be 
summarized  with  reference  to  Fig.  10,  11,  12,  that  show  diagrams  obtained 
with  the  data  given  for  the  antenna  referred  to  by  BLPZZ  19c  in  Greiser's 
and  Maye's  work.   This  antenna  operates  from  500  to  2500  MHZ.   Fig.  10, 
for  2000  Mcs  is  what  can  be  called  the  typical  pattern.   Most  computer 
runs  show  that  kind  of  diagram.   Half  power  beam  widths  of  the  order  of  90° 
in  azimuth  and  35°  in  elevation.   The  directivity,  computed  by  the 
approximate  formula  given  by  Kraus  [6],  is  around  7.5  db.   Also,  in  almost 
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Fig.    10   -   Radiation   pattern 
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Fig.  11  -  Radiation  pattern 
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Fig.    12  -  Radiation  pattern 
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every  case  the  difference  between  the  results  for  the  two  types  of 
current  distribution  is  small,  although  the  distribution  of  curve  1  in 
Fig,  9  gives  slightly  better  results,  particularly  in  what  concerns  the 
front-to-back  ratio. 

After  an  extensive  set  of  runs  was  obtained,  it  was  found  that  for 
certain  frequencies  the  diagrams  depart  appreciably  from  what  was 
expected.   Two  extreme  examples  are  shown  in  Fig.  11  and  12.   The  first 
is  "too  good,"  and  the  second  is  "too  bad."  The  most  likely  explanation 
that  we  can  find  for  that  is  in  the  way  the  current  distribution  is 
processed  in  the  program.   The  program  has  a  sample  of  the  curves  in 
Fig,  9  and  assigns  a  value  to  each  element.   Now,  that  sample  is  always 
the  same,  but  the  real  situation  in  the  antenna  varies  with  the  frequency. 
When  it  happens  that  the  wires  forming  one  element  have  length  of  one 
quarter  of  a  wavelength  the  sample  fits  exactly,  but  when  that  dimension 
"falls"  somewhere  between  two  elements,  it  is  clear  that  we  are  doing  a 
rough  approximation,  because  in  that  case  the  program  takes  the  element 
just  below  it,  and  treats  that  element  as  if  it  were  a  quarter  of  a  wave- 
length. 

To  explore  this  point  further,  we  made  a  set  of  runs  at  closely 
spaced  frequencies,  covering  one  period  of  frequency  values.   Table  I 
shows  the  results.   For  each  frequency,  column  2  shows  the  corresponding 
quarter  of  a  wavelength,  columns  3  and  4  show  the  sizes  of  the  wires 
nearest  to  that  value,  columns  5  and  7  show  the  obtained  directivity  and 
orientation  of  the  E  vector  all  around  the  ground  plane  for  the  current 
distribution  of  curve  1  in  Fig,  9,  and  columns  6  and  8  show  the  same 
results  for  the  current  given  by  curve  2. 
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1 

2 

3 

4 

5 

6      ! 

7 

8 

Freq 

A/4 

Wire1 s 

length 

Current  Dist 

Polarization 

(1) 

(2) 

MHz 

DIR/ISO 

DIR/ISO5 

(1)         (2) 

820 

.091 

.078 

.092 

5.24 

5.00 

-90° 

I 

-90° 

855 

.088 

.078 

.092 

1.45 

7.05 

+90° 

+90° 

875 

.086 

.078 

.092 

6.36 

6.06 

+90° 

+90° 

895 

.084 

.078 

.092 

5.24 

3.9 

+90° 

+90° 

935 

.080 

.078 

.092 

8.89 

10.9 

+90°within 
+60° 

+90° remain, 
directions 

-90°  within 
+60° 

+90°remain. 
directions 

970 

.077 

.067 

.078 

7.05 

7.45 

-90°within 
+135° 

+90° remain, 
directions 

-90°within 
+130° 

+90° remain, 
directions 

Table  I 
Variations  in  Results  Along  One  Period  of  Frequency  Values 
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Here  it  becomes  clearer  that  the  current  distribution  1  results  in 
a  more  balanced  set  of  results  than  distribution  2.   It  is  interesting 
to  see  the  way  the  direction  of  the  field  vector  changes  along  the  period. 
Actually,  this  change  of  frequency  is  equivalent  to  moving  the  structure 
one  period  in  space,  at  fixed  frequency,  and,  in  any  case,  we  would 
expect  the  change  in  phase  that  has  been  observed  in  all  LP's.   It  has 
been  called  the  phase  rotation  principle — the  phase  of  the  resultant  field 
changes  by  360^  with  a  change  of  one  period,  in  space  or  frequency. 
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2.3   The  BLPZZ  Antenna  as  a  Non-Uniform  Transmission  Line 

With  the  purpose  of  obtaining  a  better  description  of  the  current 
distribution  in  the  BLPZZ,  we  consider  the  structure  and  its  image  as 
forming  a  non-uniform  transmission  line.   For  easier  handling  of  the 
equations  involved,  we  split  the  structure  in  sections,  each  composed  of 
a  radiator  element  and  the  following  connecting  stub.   Then  we  apply  to 
the  radiator  Schelkunof f ' s  equations  for  slightly  non-uniform  lines, 
and  assume  that  the  stubs  are  uniform  lines  with  characteristic  impedance 
equal  to  their  average  impedance.   To  account  for  radiation,  we  introduce 
in  the  equations  an  attenuation  factor  and  a  phase  factor  obtained  from 
the  radiation  resistance  of  each  element  and  its  average  impedance.   The 
value  of  radiation  resistance  used  is  one  obtained  by  Chaney  [24]  and  Klamm 
[25],  for  rhombic  antennas  in  free  space. 

According  to  Schelkunof f  [2]  a  non-uniform  transmission  line  can  be 
described  by  equations  of  the  form 


V(x)  =  VQ(x)  +  V1(x)  +  V2(x)  +  ...  (25a) 


I(x)  =  IQ(x)  +  I1(x)  +  I2(x)  +  ...  (25b) 

where  Vn(x)  and  In(x)  are  the  expression  for  a  uniform  line  with  the 
average  characteristics  of  the  line  under  consideration.   When  a  line  or 
a  section  of  line  is  slightly  non-uniform  the  terms  of  order  2  and  above 
can  be  neglected,  without  large  error.   In  particular,  when  the  line  has 
the  property  of  having  a  constant  ZY  product,  where  Z  is  the  line 
reactance  and  Y  is  the  line  susceptance,  and  Z  and  Y  vary  slowly  about 
their  average  values,  the  first  order  term  is  given  by 
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V,l*)  rV„LBM  COS^Px   -/\(x)smhrx]  -  K„  I0[/U*)CoshPx-B(*)  Sir. A  Px]         (26a) 


I,  Wr-^J  [B(«)siT./»rx-AWCosliP*]   t   I0[AW  Sir^PK-BOOcoshRf]       (.^j 
Ko 


where 

T  =  /ZY    is  the  propagation  factor 


K  =/Z  /Y   is  the  average  impedance  of  the  line 


1    x  * 
A(x)  ■  ±-  J   Z  (s)  cosh  2rs  ds 

K0  0 


B(x)  -  ~-  J   Z  (s)  sinh  2Ts  ds 


x  A 

f 
^0  0 

z  =  z  -  z 

av 

ZA  =  Z     is  the  average  reactance  of  the  line 
0    av  6 

Y_  =  Y     is  the  average  susceptance  of  the  line, 
and  we  assume  that  the  wires  are  perfect  conductors,  perfectly  isolated, 
so  that  there  is  no  resistance  or  conductance  to  account  for. 

The  required  condition  of  slow  variation  of  line  parameters  can  be 
met  by  considering  conveniently  small  sections  of  line  at  a  time. 

Using  the  hyperbolic  form  of  the  equations  for  a  uniform  line,  we  have 

V  (x)  =  V  cosh  Tx  -  KQI0  sinh  Tx  (27a) 

V0 
In(x)  -  -  -r-  sinh  Tx  +  I  cosh  Tx  (27b) 

0        KQ  0 
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and,  then, 


V(xUVj(l^B()())coSM^-A(><)^^rVJ-|<oI0L/A(x)COShrx+(|-f3(X))si-nKn<] 

(2U) 


I(x)=^L(,4-B(x))sm/iPx-/\()f)cos/,Pj+I0[A(x)sir1hrx^(!-B(x))  coshPx] 
K°  (2  8  b) 

For  the  line  parameters  we  use  the  equations  for  the  inductance 
and  capacitance  per  unit  length  of  a  two  wire  line,  in  the  logarithmic 
form,  implying  the  assumption  that  the  radius  of  the  wires  is  very  small 
compared  with  their  spacing,  and,  so, 

Z  =  j  2H  In  (  -  )  =  j  120k  In  (  -  )  (29) 

J  7T       a  a 

Y  =  j  w  e  u/ln  (  -  )  =  j  k/120  In  (  -  )  (30) 

Si  3. 

where 

w  =  2tt  times  the  frequency 

u  =  4tt  x  10   Henry/m  is  permeability  of  free  space 

-9 
e  =  10   /36tt  Farad/m  is  permittivity  of  free  space 

s  is  the  spacing  between  wires,  and 

a  is  the  wire  radius. 
The  propagation  factor,  T,  is  a  constant  if  we  do  not  account  for 
losses.   If  we  consider  the  radiation  resistance  as  a  distributed  loss 
resistance,  with  value  Ra  per  unit  length,  we  can  use  the  following 
approximations  described  in  the  literature  for  lossy  transmission  lines, 
and  define  an  attenuation  per  unit  length. 

a  =  -Z7T-  (31) 

and  an  average  phase  factor 
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where 


8Q   '       1     8Q  / 


R 

zo 


k  is  the  free  space  phase  constant,  and 
r  =  a  +  j3 

As  mentioned  before,  we  use  for  the  computation  of  radiation  resist- 
ance the  equations  given  by  Chaney  and  Klamm  for  a  rhombic  antenna  in 
free  space.   Those  equations  were  developed  for  an  antenna  well  terminated, 
and  carrying  a  travelling  wave.   Actually,  the  present  structure  does  not 
exactly  have  the  rhombic  shape,   but   the  elevation  angles  considered  are 
small  enough  to  allow  the  assumption  that  they  can  be  neglected.   More 
daring  is  the  assumption  of  a  travelling  wave.   However,  from  the  work 
reported  by  DuHamel  on  nonuniform  transmission  lines  [16],  we  conclude 
that  the  value  of  Q  should  not  be  expected  to  be  critical,  and  that  is 
essentially  what  would  change  directly  with  the  value  of  radiation  resist- 
ance. 

So  far  we  have  not  accounted  for  the  mutual  impedance  among  radiators, 
and,  naturally,  that  is  due  to  the  large  amount  of  computations  needed, 
that  would  exceed  the  time  allowed  for  the  present  work.   After  presenting 
the  results  obtained  with  the  present  assumptions  we  will  try  to  do  a 
rough  estimate  of  the  mutuals  effects,  based  on  the  values  found  for 
current  distribution  and  impedances. 

Using  equations  (28)  for  the  radiators  and  equations  (27)  for  the 
stubs,  we  start  from  the  end  of  the  line,  where  there  is  a  short-circuit, 
and  proceed  toward  the  feed  point  computing  the  input  impedances,  and  the 
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coefficients  of  the  equations.   As  a  block  the  equations  can  be  expressed 
by 


V(x)  =  A  VQ  -  KQ  B  IQ  (33a) 


I(x)  =  -r0vo  +  Dlo  (33b) 

and,  for  each  section,  we  obtain  A,  B,  C  and  D,  and  the  ratio  V  /I  . 
When  we  reach  the  feed  end  of  the  line  ratio  V-/I-  is  its  input  impedance. 
Then  we  assign  to  the  current  the  value  1.0  +  jO.O  and  we  proceed  forward, 
using  the  previously  computed  values  of  A,  B,  C  and  D,  obtaining  the 
values  of  current  and  voltage  along  the  line,  for  as  many  points  as  the 
number  of  sections  considered.   From  these,  and  using  the  equations,  the 
complete  description  of  the  line  can  be  obtained. 

We  program  these  equations  into  the  IBM  360  system,  according  to  the 
program  in  Appendix  C.   Appendix  D  shows  the  details  of  the  equations 
involved. 

Sample  program  SMIGUEL  shows  the  results  for  the  same  antenna  that 
was  used  in  sample  program  AZORES  (radiation  pattern) ,  for  three  frequencies 
within  the  operation  band — 1100,  1700  and  2300  MHZ.   For  each  frequency 
the  first  run  does  not  account  for  mutual  impedance  effects.   The  second 
run  includes  the  mutual  resistance  that  would  be  found  between  parallel 
wires  with  sinusoidal  current  distribution. 

The  experimental  results  obtained  by  Greiser  and  Mayes  [8],  [9],  with 
the  BLPZZ  fed  against  ground  show  that  the  input  impedance  has  a  resistive 
component  varying  between  90  and  250  ohms  and  a  reactive  component 
varying  between  +10  and  -100  ohms.   Since  we  are  considering  a  double 
structure,  we  would  expect  the  resistance  values  to  be  about  doubled. 
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Although  the  behavior  of  the  reactive  component  is  more  difficult  to 
predict,  we  would  expect  the  same  orders  of  magnitude. 

These  assertions  are  confirmed  in  the  results  obtained  for  2300  MHz. 
Referring  to  SMIGUEL  the  column  of  ZINR  values  shows  the  input  impedance 
seen  at  the  feed  end  of  each  radiator,  and,  not  only  the  input  impedance 
is  as  expected  (291. 3-j314,8  ohms),  but  also,  as  assumed  before,  each 
element  essentially  provides  a  match  for  the  next  element.   The  column  of 
CUR  values  shows  the  magnitude  of  the  current  at  the  feed  end,  at  the 
vertex,  and  at  the  stub  end  of  each  radiator.   We  note  that,  for  this 
frequency,  element  2  is  smaller  than  a  quarter  of  a  wavelength,  and  that 
element  3  is  larger  than  that.   Up  to  the  third  element  the  current 
magnitude  values  are  typical  of  an  attenuated  standing  wave  situation. 
After  element  3  the  distribution  is  of  a  heavily  attenuated  travelling 
wave. 

As  mentioned  beofre,  the  more  delicate  assumption  in  this  approach 
was  to  apply  to  this  structure  the  radiation  resistance  given  by  Chaney 
and  Klamm  [24],  [25]  for  the  rhombic  antenna.   It  appears  from  these 
first  results  that  such  assumption  is  reasonable  for  the  large  elements, 
but  it  is  not  adequate  for  the  small  elements.   Also,  for  these  elements, 
the  effects  of  mutual  impedance  become  increasingly  important,  since 
they  are  "electrically  closer"  to  each  other.   These  two  sources  of  error 
become  more  important  as  we  lower  the  frequency.   At  1700  MHz  (the 
resonance  dimension  falling  between  elements  7  and  8)  the  results 
obtained  for  the  first  4  elements  are  inconsistent. 

A  closer  inspection  of  the  results  shows  that  the  departures  in 
impedance  values  are  related  to  extremely  high  values  of  Q  (above  20),  and 
that  fact  even  more  strengthens  the  idea  that  the  loading  due  to  mutual 
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effects  plays  an  important  part  in  the  behavior  of  the  "feed  region"  of 
the  antenna. 

At  this  point  we  found  it  useful  to  make  some  numerical  analysis  of 
radiation  and  mutual  resistance  values,  with  the  following  conclusions: 

(1)  The  radiation  resistance  of  the  "rhombic"  elements  increases 
steadily  with  size,  from  a  few  ohms  at  the  "feed  region"  to 
thousands  of  ohms  at  the  "end  region," 

(2)  The  mutual  resistance  between  parallel  wires  in  adjacent 
elements  with  sinusoidal  distribution  of  current,  in  phase,  decreases 
steadily  from  hundreds  of  ohms  at  the  "feed  region"  to  essentially 
zero  at  the  "end  region." 

This  enables  us  to  obtain  a  better  approximation,  by  also  including 
the  mutual  effects,  taking  two  parallel  consecutive  wires  of  the  structure 
at  a  time.   Again  here,  we  neglect  a  few  points;  namely,  that^  the  parallel 
wires  considered  do  not  have  the  same  size,  the  input  currents  are  not  in 
phase,  the  current  distribution  is  not  exactly  sinusoidal,  and  that  there 
are  other  interactions  to  be  accounted  for,  among  the  non-parallel  wires. 
However,  before  taking  this  step  we  made  extensive  numerical  calculations 
with  different  arrangements  of  wires  with  different  conditions  of  phase 
differences  and  current  distributions,  and  the  results  obtained  were  all 
within  the  same  orders  of  magnitude. 

Essentially,  what  we  are  doing  here  is  to  introduce  a  rough  correction 
to  lower  the  Q  of  the  "feed  region,"  and  the  same  result  could  be  obtained 
if  we  assigned  to  it  specific  (and  reasonable)  values.   As  shown  in  the 
referenced  paper  by  DuHamel  [16]  the  exact  value  of  Q  is  relatively 
unimportant. 
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We  note  that,  by  this  procedure 

(1)  at  the  "feed  region"  we  are  essentially  neglecting  the 
travelling  wave  approach  since  the  rhombic  radiation  resistance  is 
much  smaller  than  the  mutual  resistance,  and 

(2)  inversely,  at  the  "end  region,"  we  are  neglecting  the  sinus- 
oidal distribution  approach. 

The  results  obtained  with  this  procedure  apparently  prove   its 
soundness.   For  the  three  cases  previously  considered,  the  input 
impedance  values  now  fall  in  the  expected  range  (439. 1-J148. 5,  562.1-J109.8, 
and  516.4+jl9,9  ohms).   The  current  distribution  shows  the  same  features  as 
before  but  with  stronger  attentuation,  making  it  resemble  more  closely  an 
attenuated  travelling  wave  configuration. 

The  next  step  would  be  to  obtain  a  more  exact  description  of  the  line 
by  considering  a  larger  number  of  smaller  sections,  obtaining  better 
approximations  for  the  non-uniformity  corrections,  and  a  more  detailed 
picture  of  the  current  values.   That  would  enable  the  application  of  these 
results  to  the  radiation  equations.   Unfortunately  the  attempts  to  program 
this  detailed  analysis  failed  due  to  unknown  reasons,  related  to  the 
difficulties  encountered  during  the  "adaptation"  phase  of  the  new  IBM  360 
system,  recently  installed  at  the  Naval  Postgraduate  School.   Although  we 
have  received  extensive  help  from  the  Computer  Facility  staff,  and 
undoubtedly  the  program  would  eventaully  work,  there  was  not  time  available 
to  pursue  the  project. 

Before  closing  this  subject  we  note  that,  due  to  the  method  used  in 
deriving  the  coefficients  in  equations  (33),  there  is  a  cumulative  source 
of  error  in  the  approximations  used,   This  is  so  because  every  new  value 
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is  based  on  another  value  that  had  been  obtained  by  an  approximation, 
and  only  the  first  value  is  an  exact  value.   Since  we  start  from  the  end 
of  the  structure,  the  "end  region"  coefficients  are  little  affected,  but 
the  ones  at  the  "feed  region"  are  likely  to  be  farther  away  from  exactness. 
This  might  be  another  reason  for  the  worse  results  obtained  with  the 
smaller  elements.   The  way  to  correct,  or  at  least  to  minimize  this  error, 
would  be  to  split  the  line  in  more  sections,  as  was  attempted. 

Finally,  we  note  that  the  current  distribution  does  not  show  the 
marked  peak  near  the  resonance  region  obtained  with  parallel  feed  type  of 
structures,  or  circuits.   On  one  side  this  seems  a  reasonable  result  due 
to  the  different  physical  lay-out  of  the  antenna;  but,  on  the  other  side, 
it  remains  to  be  explained  why  the  radiation  patterns  show  acceptable 
results  with  different  current  distributions ,   An  inspection  of  the  phase 
distribution  of  the  current  (column  CANG)  shows  that  in  all  cases  the 
current  at  the  feed  point  of  any  two  consecutive  elements  is  very  nearly 
in  phase  opposition,  except  for  the  element  with  dimensions  closest  to  the 
resonance  dimension,  and  the  one  before  that.   For  these,  the  phase 
difference  is  near  90°,  as  also  happens  in  other  LP  structures-   This  means 
then  that  only  the  "active  region"  is  important  for  radiation  characteristics, 
and  the  two  models  of  current  distribution  previously  used  are  essentially 
equivalent  to  each  other  and  to  the  computed  values. 
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APPENDIX  A 


I.   Geometric  Relations  in  the  BLPZZ 

With  the  aid  of  Fig-,  A-l  the  following  relations  are  evident: 

Ll  bl 

L  =  — J  and     b   =  — ±r  (A-l) 

n    n-i  n    n-i 

P  P 

The  relation  between  the  length  of  the  radiating  wires  and  the  stubs 
can  be  obtained  in  the  following  way: 

2b  cosy  2b1   cosy 

L  =  -7s —         L  '  =  — — (A-2) 

n    tan  aE  n     tan  as 


2b   cosy   2b   _  cosy 

L  -L   ,  =  -*■ —  - =  2b   ..siny  +  2b  siny  +  4b   ..  'siny 

n  n-1   ■  tan  aE     tan  aE        n-1         n         n-1 


cosy  b  . 

- —  (  —  -  1)  =  b   ..  siny  (1  H )  +  2b    '  siny 

tan  aE     p         n-1     '      p       n-1 


I^±S  (l-p)  -  bn.1  tanY  (14,)  +  2b   ■  P  tanY 


n    2p  tan^y  tan  aE 


b  '    (l-p)  -  (1+p)  tany  tan  aE 


(A-3) 


Other  useful  relations  are 


2L'    =  L   ,  +  L   ,  tan  aE  tany  +  L  -L  tan  aE  tany 
n-1    n-1    n-1  n  n 


L  ,  (1  +  —  )  +  L   ,  (1  -  —  )  tan  aE  tany 
n-1      p      n-1      p 


L 

-^-^  [(1+p)  +  (p-1)  tan  aE  tany] 


56 


a)  Basis  structure 


b)  Bent  structure 


Fig  A-l  -  Geometric  Relations  in  the  BLPZZ 
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from  which 

L 


[(1+p)  -  (1-p)  tan  aE  tany]  tan  as 


n-1   4p  cosy 

and  combining  with  (A-3) 

tyvj _        2£>  tartf  tarnxe. 

2b      .,    cosy 

,         .  T  n-1 

and,    since  L      .,    = 

n-1  tan  aE 

{\-\>)  -(n-p)tar)$  tano(e    -  {Up)  tantf  tar) o<s  -(hf>)  tan  <*e  tan  o<stanz^ 


(  ±t-£- )  (fono<£  +  tancx5)        I  (— rj|)  (tan«(e  +tar> o<sf-4  tar)#£  tar)<xs 
Ur>t;_lzt :         ' (A_4) 

2    '  """Yf-  I4r)0<s 


Since  the  structure  is  usually  defined  by  the  parameters  aE,  as,  ,, 
b  for  the  last  element  and  the  number  of  elements,  equation  (A-4)  gives 
the  next  parameter,  y,  and  equations  (A-l)  and  (A-2)  define  the  rest  of 
the  structure.   The  length  of  the  stubs  can  be  obtained  with  (A-3),  or 
approximately  by 

b 


n  tan  as  , .  cN 

b   =  — (A-5) 

n    r~     tan  aE 


This  equation  is  exact  when  as  =  aE,  but  it  can  be  used  as  a  good 
approximation  for  the  typical  values  used.   The  equation  is  obtained  by 
inspection  of  Fig.  A-l,  with  the  help  of  a  sort  of  "optical  illusion." 
In  fact,  when  as  ■  aE  the  structure  repeats  itself  at  /p  increments,  with 
an  upside  down  rotation.   When  as  ^  aE,  however,  the  rate  of  change  of  b' 
with  as  is  not  linear  and  the  equation  is  not  true. 

The  relations  concerning  the  "phase  center"  of  the  element  are 
obtained  from  Fig.  3,  by  inspection 
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tan  &,  -     L*>  s'-n6+  b>  cos  (if- e)  (A_6) 

6^  si*  (X-l) 


tar,  S2  _    Ly,  sine  *  i>nCosC^+e) 


Clm     kj^fl)  (A.8) 

Cos  61 


C.     6_^pi  (A.9) 


II.   Transformations  of  Coordinates 

Expressions  for  the  angle  between  a  wire  and  a  direction  -9'. 
Considering  Fig.  A-2,  that  is  a  planif ication  of  the  Y=0  and  Z=0  planes, 
we  obtain  the  following  relations: 

cos  a  sin  b  =  cos6 

cos  a  cos  b  =  sin6  cos  0-        0..  =  180-0 

R  cos6'  =  R  cos  a  cos  (8-b) 

=  R  [cos  3  cos  a  cos  b  +  sin  6  cos  a  sin  b] 
=  R  [cos  8  sin  0  cos  0..  +  sin  6  cos  8] 
=  R  [sin  3  cos  0  -  cos  3  sin  0  cos  0] 
The  wire  in  the  figure  corresponds  to  wire  1  in  a  BLPZZ  element,  and 
the  angle  3  is  given  by  (see  Fig.  A-l) 

3  -  e  +  90°  -  y  =  90°  -  (y-e) 
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Fig*  A-2  -  Transformation  of  coordinates 
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so   that,    for  wire   1 

Cos6,'  =    cos  (y-£)  Co5&   -sin  (ft- a)  Sin  6  cos  0 

(.A— lv) 

For  wire  2  we  have  instead  of  an  angle  3,  and  angle  -6,  given  by 

-6  =  -  (90  -  (y+e)) 
and  using  the  equation  for  cos  6'  we  get 

cos  6>£  _  _  cos  (X*e)  cos  6  -  sir,  ($+£)  sm  $  cos  jt>  (A-ll) 

For  wire  3,  we  obtain  (see  Fig.  A-3) 
cos  a  sin  b  =  cos  0 
cos  a  cos  b  =  sin  6  cos  0 

cos6'  =  cos  a  cos  ((3-b) 

=  cos  a  cos  b  cos  6  +  cos  a  sin  b  sin  3 
=  cos  9  sin  6  +  sin  6  cos  3  cos  0.. 

and,  with  0=0     and     3  =  90  -  (y-c) 

cos6'  =  cos  0  cos  (y-e)  +  sin  0  sin  (y-e)  cos  0       (A-12) 

For  wire  4,  we  note  that  the  only  difference  is  again  the  direction 
of  the  angle  3  =  90  -  (y+e),  so  that 


cosO!  =  -cos  0  cos  (y+e)  +  sin  0  sin  (y+e)  cos  0      (A-13) 
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Fig,  A-3  -  Transformation  of  Coordinates  -  Wire  3 
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III.   The  Components  of  the  E  '  Vector  in  Spherical  Coordinates 

According  to  Fig.  A-4,  we  obtain  the  following  relations,  for  wire  1 

E0,  =  fe,'  C05/4'  (A-14a) 

E0,    =  £9l  5'nA'  (A-14b) 

Using   the   expression  from  spherical   trigonometry 

,  cos   a  -  cos   b   cos   c 

cos   An    =  ; : : 

1  sm  b   sin  e 

we  obtain 


cos  A   -    cos  (4o- p)- cos  $',  cose        siT)p-co&d(smficos$-cospsin0cos0) 

sm  e;  sm  e  sjhcofe;    sm  a 


sin  ft-  smftcos'e  t  cosp  Sin0cos$cosf6       S/n/3(hcos  6)+Cosp>S'm&cosScosfi 
\j  I- cos1 9',     $'i*6  \j\-cosl9',    sirtB 


sm  p  sin  $  f  cosj3  Cos  $  costf 
]/  !•  cosiO'l 


and,    since   3   =   90  -    (y-e), 

cos  A   _     cos  Or- a)  sm  9  +  sin  (f-£)  cos$cos  tf 


I (A-15a) 

\/  i-  cos2e] 


We  have  the  denomination  in  the  form  \l  l~Co$  #/  because  in  the  field 
expression  we  have  a  sin  8'  in  numerator  that  cancels  with  this.   Now  we 
obtain,  from  here,  sin  A 
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Fig.  A-4  -  Components  of  the  field  vector 
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SmzA   =  I- cos1  A    _     f-(smZfiCosi6+cos*j&siril6cos*0-2Sir)pcospSmdcosdcosfi) 


(S/nZj&S/r)*&+Cos*£  Cosz0  cesl0  +  Z  $mf$cos>p  S>n6 Cosd cOs0) 


that  reduces  to 

StrfA,  .     j-Sivfp-tos^O-Smp)  _    (\-$m*p)(\-Co<?0)    _       cosl/SSm'lfi 
Sin*&l  sin*  0',  Sin*  B\ 

and 

Sir,  A,   -      cosptinj    ,       $ir>(X-S)sin&  (A-15b) 

Srn  0[  sh  £>{ 

The  procedure  is  analogous  for  the  other  wires,  and  the  results  are 
for  wire  2 

E.   =  -  EQ,  cos  A.  (A-16a) 


EQ  =  Ee,  sin  A2  (A-16b) 


COsAi    .     Cos  (X+E)  sin  9-  Smti  +  8)cos9cos<f 

$ir>  $1  (A-17a) 


Sir,  A2    -      sin  (if+e )  Sim  0  (A-17b) 
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for  wire   3 


EQ      =   Ee,    cos   A3  (A-18a) 


E.     -  -  E   ,    sin  A  (A-18b) 

W3  3 


CosA^  -     Co s  (f-e)  stn  B-  sin  (ITS)  cos  &  cos0  (A-19a) 

Sir,  ^ 

Sir)A3   -  Srn(r-e)  $>"  0  (A-19b) 

Sin  0^ 


and  for  wire  4 


E.   -  -  Eof  cos  A.  (A-20a) 

64      64      4 


E0  =  E0,  sin  A4  (A-20b) 

4     4 


cos/4*-  cos(t+e)sm*  *sin(!+e)cos&cos0 

sin  /)4  =  si*(ir+e)sm0  (A-21b) 

IV.   The  Angle  a  for  the  Correction  of  Reference  Point 

This  is  the  angle  between  the  direction  of  c  (see  Fig.  3)  and  the 
direction  to  a  point  in  the  far  field.   Since  c  is  in  the  same  vertical 
plane  as  the  wires,  we  can  use  the  relations  developed  for  the  conversion 
of  the  angle  0'  to  spherical  coordinates,  with  6  in  place  of  6,  and  we 
obtain  the  following  relations: 

cos  a,  =  -[sin  6..  cos  6  +  cos  6.  sin  9  cos  0J  (A-22) 

cos  a„  =   sin  61  cos  6  -  cos  5.  sin  9  cos  0  (A-23) 

cos  a~  =   sin  5.  cos  0  -  cos  S..  sin  9  cos  0  (A-24) 

cos  a,  =   sin  6_  cos  9  +  cos  6~  sin  9  cos  0  (A-25) 
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Fig.  A-5  -  Reference  point  and  direction  angles 
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The  angles  6±   and  52  are  shown  in  Fig.  3,  and  given  by  equations 
(A-6)  and  (A-7).   Fig.  A- 5  shows  the  way  how  the  a   angles  are  defined. 
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APPENDIX  B 

Field  Equations  for  Wires  2,  3  and  4 
I  -  Wire  2 

Since  wire  2  (see  Fig.  5)  is  in  series  with  wire  1,  we  need  to  add  in 
the  exponent  a  factor  corresponding  to  the  length  of  wire  1  -  2b  -  and  we 
get  from  (7) 


E„,  „  if  t  e  _£ =^e_  (B-D 


Now,  using  equations  (A-ll) ,  (A-16),  (A-17)  and  (A-23) ,  we  obtain 


EB  --it  1,  e         K  _e -e „  ,„  ,, 

e>"    *  [-tik(^. .«..«)]  <B-2> 

s/y)  ($*-£)  sin  tf 
correction  fo: 

being  c9   and  cos   ou   given  by   equations    (A-9)    and    (A-23). 


~ ■  "1  k  C    COS  GL 

where  the  factor  e     2      2   is  the  correction  for  the  reference  point, 


II  -  Wire  3 

Here  we  consider  the  current  flowing  upward,  but  to  use  the  same 
exponential  representation  for  the  current,  we  consider  the  distance  as 
negative  away  from  the  source.   This  is  the  usual  representation  for  trans- 
mission lines,  inherent  in  the  present  assumption  of  a  travelling  wave  in 
the  antenna.   Then 
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I  =  I  e"j6s 
3    0 


and 


E  ,  Fk  I,  sft,  4  e    e -e     (b-4) 

Now,  using  equations  (A-12) ,  (A-18) ,  (A-19)  and  (A-24),  we  obtain 


_fr[f-t>+<±c**«*]   HHIf^)]b  -NKf^)]b 


x  [<^5  (/-£)  5/Vj  d  -  Sf'n  (#-6)  cos  &cos0] 


X 
(B-5) 


^[p+c3cos«,]     fa^+ostyL    -[**iM-f«**0> 

>< 
(B- 


t^--^]»e  _e . ^_£ x 


^   Sin  (%-t)  i'n  0 


wh 


ere  c„  and  a„  are  given  by  equations  (A-8)  and  (A-24).   Note  that 


C3  =  CV    and  °4  =  °2  ^Fig#  5^* 


III.   Wire  4 


This  wire  is  in  series  with  wire  3,  so  that  the  procedure  is  analogous 
to  wire  2  that  was  in  series  with  wire  1 


fce' "~  z  ^  sm  °i e         *-  ■  -  ~"e  (B_7) 


Now,    using  equations    (A-13) ,    (A-20) ,    (A-21)    and    (A-25),   we  obtain 


X  \cos(p£.)  Sm8  +  5m  (Jff£j  ioj(9  coi^j 

(B-8) 
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(B-9) 


5/n  (J+£)  5/rj  ^ 
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APPENDIX  C 

Equations  for  Non-Uniform  Line 
I.   Line  Parameters 

Neglecting  wire  resistance  and  conductance  between  wires,  the  line 
parameters  are 

Z  =  j  120  k  In  (s/a) 
and 

Y  =  j  k/  (120  In  (s/a)) 
as  given  by  equations. 

From  here  we  obtain  the  line  impedance 

K  =  Jzjx   =  120  In  (s/a)  =  fr 

and  the  propagation  factor  not  accounting  for  losses 

r  -  /ZY  =  jk 

In  terms  of  the  structure  defining  dimensions  and  angles  (see  Fig. 
A-la) ,  we  obtain  the  following  equations 
For  the  radiators 

Z(x)  =  j  120k  (B  +  In  (A  +  x) )     for  0  <  x  <  2b     (C-la) 
Z(x)  =  j  120k  (C  +  In  (D  -  x))     for  2b  <_  x  <  4b    (C-lb) 


where 


1  -  2b  siny 

A  =  7 s —  sine 

cos  (y-e; 


B  =  ln  2  cos  (y-s) 
a 
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C  «  in  l_££s_Mel 


L  -  2b  sin  y  gine+  2b__cos_Jxzil  +  2b 
cos  (Y+e)  cos  (Y+O 


For  the  stubs 

z(x)  =  j120kT  E  +  In  (P+x)  J 


where 


E  =  In 


2sin  Ysin  e 


2b  tan  or 


F  = 


s 


7p~sin  y     yp~tan  a£ 
and  the  quantities  involved  are  as  defined  in  Fig.  A-la 


(C-2) 


*.n<?,i)  :*,S^) 


2L.M 


'■$  (r,vt-»)      *J&)        zjtfj) 


Fig.  C-l    Input  Impedance, 
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II.   Radiation  Resistance 

According  to  Chaney  and  Klamm  [24],  [25],  the  radiation  resistance  of  a 
rhombic  antenna  in  free  space  is  given  by 

Rr  =  120  |l.l544  +  2  In  (2  kl  sin2  a)  +  2  Ci  (2kl)  -  2  Ci(2kl  sin  a) 

-  Ci  (2kl  (1  +  cos  a)  -  Ci  (^2kl  (1  -  cos  a)  J 


f    2kl    sin   aj     Ci  (2kl    cos   a    (1    -   cos   a)    J  +   Ci  Clkl    sin   a    (1    -    sin   a)    J 


+   cos  (   2kl  sin  a 


+  Ci  (2kl  cos  a(l+cos  a)  J  +  CiT2  kl  sin  a(l+sin  a) J  -  2Ci  f   2  kl  cos   a  J 


C   2kl  sin2  a  ^ 


2Ci   2kl  si 


+  sin  (   2kl  sin  a)        Si  (   2kl  cos  a  (1-cos  a) J  -    Si  f  2kl  sin  a  (1-sin  a)  ) 


Si  ^2kl  cos  a  (1+cos  a) J  +   Si  (^2kl  sin  a  (1+sin  a)  J  +  2  Si  (  2kl  cos   a  J 


(^  2kl  sin2  a  ")   j 


-  2  Si   2kl  si 


where  1  is  the  wire  length  -  2b  -  ,  a  is  the  vertex  angle  90   -  y 
Ci  is  the  cosine  integral 


ci(x)  = 


COS  V 


dv 
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and  Si  is  the  sine  integral 

x 
Si(x)  =  /  ^B^L   dv 

0    v 

III.   Input  Impedance  and  Current  Distribution 

Using  equations  (33)  and  dividing  the  line  in  sections  as  shown  in 
Fig.  D-l,  we  obtain  for  each  section 


Z.  (I,n)  =  Kave(I.n)  ^'"K^^  V^1?^  iff'" 
in  Zxn(I,n+l)  C(I,n)  +  Kave(I,n)  A(I,n 


Each  stub  counts  as  only  one  section  of  uniform  line.   The  number  M 
of  divisions  within  each  radiation  depends  on  the  degree  of  approximation 
desired,  but  it  is  convenient  to  use  for  M  an  even  number,  to  avoid  the 
relatively  faster  rate  of  change  of  characteristics  at  the  middle  point. 
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APPENDIX  D 

Computer  Programming 
I.   Program  AZORES 

This  program  is  written  in  code  FORTRAN  IV  and  it  was  prepared  for 
the  IBM  360  system  that  is  operating  at  the  Naval  Postgraduate  School. 

The  program  uses  two  subroutines.  The  main  program  reads  the  data, 
computes  and  writes  the  physical  characteristics  of  the  antenna,  selects 
the  current  distribution  for  the  antenna  wires,  controls  the  coordinates 
of  the  desired  pattern,  prints  the  results,  and  calls  for  the  subroutines. 

The  pattern  computation  is  made  by  subroutine  CORVO,  that  is  given 
the  antenna  parameters,  the  elevation  and  azimuth  angles,  and  the  current 
distribution.   In  output,  CORVO  gives  the  two  components  of  the  electric 
field. 

The  program  computes  each  pattern  twice,  for  two  different  models  of 
current  distribution.  Model  1  (see  Fig.  9)  is  fed  in  data  cards;  Model  2 
is  generated  by  the  main  program. 

The  use  of  a  large  number  of  comment  cards,  hopefully  makes  the 
program  self-explanatory.   As  much  as  possible,  the  names  given  to 
algebraic  and  geometric  symbols  have  the  same  phonetics  as  the  ones  used 
in  the  text.   The  only  remarkable  exception  is  the  symbol  TAU,  that  stands 
for  p,  the  structure  scaling  factor. 

The  following  pages  contain  the  programming  statements  and  the  output 
print  of  a  sample  program. 
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-15. 

636 

0  .  c 

85. 

0. 

,120 

90, 

,000 

-14. 

907 

0.0 

90. 

0. 

121 

90. 

.000 

-14. 

863 

0.0 

95. 

c 

,117 

90, 

.000 

-15. 

145 

ICO. 

c 

,  113 

90, 

.00  0 

-15. 

427 

0  .  c 

105. 

c 

,113 

90, 

.000 

-15. 

424 

CO 

110. 

c 

,  1  19 

90. 

,000 

-14. 

985 

0.0 

115. 

c 

,  131 

90, 

.000 

-14. 

153 

U.J 

120. 

c 

,  148 

90, 

,00  0 

-13. 

103 

0.0 

125. 

n 

u  « 

,167 

90. 

,000 

-12. 

0  29 

0.0 

130. 

0. 

,187 

90, 

,000 

-11. 

066 

0.0 

135. 

c 

,204 

90, 

,000 

-10. 

285 

CO 

140. 

0. 

,219 

90, 

.000 

-9. 

708 

CO 

145. 

c 

,228 

90, 

,000 

-9. 

323 

0  .  c 

150. 

c 

234 

90, 

.000 

-9. 

104 

0.0 

155. 

c 

237 

90. 

.000 

-9. 

013 

0.0 

160. 

c 

,237 

90, 

.00  0 

-9. 

Oil 

0.0 

165. 

c. 

235 

90, 

.000 

-9. 

060 

0.0 

170. 

c 

,234 

90. 

.000 

-9. 

12? 

CO 

175. 

c 

,232 

90, 

,000 

-9. 

172 

0.0 

180. 

c 

232 

90, 

,000 

-9. 

191 

LIRECT IV  I TY/ISC      FI  3CE 

i 

ELEV 

3DB 

10.969CB 

33. 

25. 

THIS  DIRECTIVITY  FIGURE  WAS  COMPUTED 

ASSUMING  THAT  THE  PATTERN  IS  WELL-BEHAVED 
IT  HAS  NO  MEANING  IF  THERE  APE  SIDE-LOBES 
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CURRENT  L 

'  T  S  TR  I  OUT  ICN  1 

FREQ  = 

0.1  I  OF  10 

hLFV 

FI 

F 

PSI 

FO 

O.C 

O.C 

3. 

,228 

90. 

,000 

-0. 

000 

1C. 

0.0 

3. 

,14  3 

90. 

,00  0 

-0. 

232 

20. 

2, 

,900 

90, 

,000 

-0. 

930 

3C. 

CO 

2, 

,537 

90. 

,00  0 

-2. 

093 

40. 

CO 

2, 

,104 

90. 

,000 

-3. 

7?0 

50. 

f  .0 

1, 

,653 

90, 

,000 

-5. 

816 

6  0  . 

c.c 

1. 

,225 

90. 

00  0 

-8. 

416 

70. 

0.0 

0, 

,84  6 

90. 

,000 

-11. 

636 

80. 

CO 

0, 

,525 

90. 

,000 

-15. 

768 

90. 

180. 

0. 

,270 

-90. 

,000 

-21. 

539 

80. 

180. 

c. 

,08  2 

-89, 

999 

-31. 

878 

70. 

180. 

0. 

,043 

89. 

999 

-37. 

435 

60. 

18  0. 

0. 

,117 

90. 

000 

-28. 

739 

50. 

180. 

c« 

,  156 

90, 

,000 

-26. 

30  8 

40. 

180. 

0. 

,174 

90. 

00  0 

-25. 

358 

30. 

180. 

0, 

,182 

90. 

,000 

-24. 

974 

20. 

180. 

c. 

,  186 

90. 

,00  0 

-24. 

793 

10. 

180. 

c. 

,188 

90, 

,000 

-24. 

690 

0.0 

18  0. 

0. 

,  189 

90. 

000 

-24. 

653 

C  .  0 

cc 

3. 

,228 

90. 

,00  0 

0. 

o 

0.0 

5. 

1 

,224 

90. 

,000 

-0. 

012 

0.0 

10. 

3, 

,211 

90. 

,000 

-0. 

^47 

r  o 

15. 

3. 

,  189 

90. 

ooo 

-0. 

106 

0  . 0 

20. 

3, 

,  15? 

90, 

,000 

-0. 

190 

0.0 

25. 

3, 

,119 

90, 

,  DOO 

-0. 

3  00 

0.0 

3  0. 

■* , 

,070 

90, 

,00  0 

-0. 

4  37 

0.0 

35. 

3, 

.011 

90. 

,00  0 

-0. 

605 

0.0 

40. 

2, 

,941 

90. 

,000 

-c. 

,8  08 

45. 

2. 

,860 

90. 

,00  0 

-1. 

,052 

0  .0 

50. 

2  . 

,  765 

90. 

,000 

-1. 

344 

0.0 

55. 

2. 

,656 

90, 

,00  0 

-1. 

694 

\j   a  ' J 

60. 

2. 

,531 

90. 

,00  0 

-2. 

114 

u  .  0 

65. 

2. 

,389 

90. 

,  000 

-2. 

616 

CO 

70. 

2. 

.230 

90, 

,00Q 

-3. 

213 

0.0 

75. 

2, 

,056 

90. 

,000 

-3. 

920 

80. 

1. 

,869 

90. 

,000 

-4. 

748 

0.0 

85. 

1. 

,673 

90. 

,  0  0  0 

-5. 

,707 

O.C 

90. 

1. 

,475 

90. 

,000 

-6. 

802 

0.0 

95. 

1. 

,280 

90. 

,000 

-8. 

0  34 

O.C 

100. 

1. 

,094 

90. 

,000 

-9. 

397 

0.0 

105. 

C 

,923 

90. 

,00  0 

-10. 

877 

0.0 

110. 

c . 

,770 

90. 

,00  0 

-12. 

449 

0.0 

115. 

n 

,63  8 

90. 

,00  0 

-14. 

079 

0  .0 

120. 

c 

,52* 

90. 

,000 

-15. 

722 

O.C 

125. 

0. 

,439 

90. 

,000 

-17, 

325 

CO 

130. 

c. 

,369 

90, 

,000 

-18. 

829 

o.o 

13  5. 

c 

,316 

90. 

-20. 

184 

140. 

0. 

.276 

90, 

,000 

-21. 

350 

0.0 

145. 

0, 

,247 

90. 

,000 

-22. 

310 

CO 

150. 

0. 

,227 

90. 

,000 

-23. 

064 

0.0 

L55. 

0, 

,212 

90. 

,000 

-23. 

633 

0.0 

160. 

c, 

,203 

90. 

,000 

-24. 

046 

CO 

165. 

0. 

,  196 

90. 

,000 

-24. 

331 

0  .  0 

170. 

0. 

,192 

90, 

,000 

-24. 

516 

CO 

175. 

c, 

,190 

90, 

,000 

-24. 

619 

O.C 

180. 

0, 

,  189 

90, 

,000 

-24. 

653 

CIRECTIVITY/ISO 
6.368CE         68. 


FI  308 


35. 


FLEV  30Q 


THIS  DIRECTIVITY  FIGURE  WAS  COMPUTED 

ASSUMING  THAT  THE  PATTERN  IS  WELL-8EHA VED . 
IT  HAS  NO  LEANING  IF  THERE  ARE  SIDE-LORES. 
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CURRINT    riSriUHUTICN    2 

FPEU  = 

O.UOE     10 

tl.LV 

FI 

E 

PSI 

E08 

0.0 

0.0 

2.951 

90.000 

CO 

10. 

0  .  c 

2.87Q 

90.000 

-0.216 

20. 

0.0 

2.675 

90.000 

-0.852 

30. 

CO 

2.377 

90.000 

-1.878 

40. 

0.0 

2.02^ 

90.000 

-3.254 

50. 

0.0 

1.664 

90.000 

-4.976 

60. 

o.c 

1.297 

90.000 

-7. 14? 

7C. 

0.0 

0.930 

9C.C00 

-10.028 

80. 

0.0 

0.580 

9C.O00 

-14. 132 

9C. 

180. 

C.284 

-90.000 

-20.345 

80. 

190. 

0.079 

-83.999 

-31.488 

70. 

180. 

0.031 

89.999 

-33.6  76 

60  . 

180. 

C.078 

90.000 

-31.517 

50. 

18  0. 

C.  109 

90.000 

-28.629 

4C 

1 8  0 . 

C.149 

90.000 

-25.949 

30. 

180. 

C.197 

90.000 

-23.494 

20. 

180. 

0.243 

90.000 

-21.678 

10. 

180. 

C.275 

90.000 

-20.611 

C.C 

180. 

0.286 

90.000 

-20.264 

0  . 0 

CO 

2.951 

90.000 

0.0 

O.C 

5. 

2.949 

90.000 

-0.006 

C  .  0 

10. 

2.343 

90.000 

-0.0  24 

J  .  0 

15. 

2.934 

90.000 

-0.052 

C  .  G 

zc. 

2.921 

90.000 

-0.038 

0.0 

25. 

2.906 

90.000 

-0. 133 

0  . 0 

30. 

2.889 

90.000 

-0.  134 

0.0 

35. 

2.869 

90.000 

-0.244 

0.0 

4  0 . 

2.845 

90.000 

-0.  H7 

0.0 

45. 

2.814 

90.000 

-0.413 

0  .  c 

50. 

2.772 

90.000 

-0. 544 

P       o 

55. 

2.713 

90.000 

-0.731 

c .  c 

60. 

2.632 

9  0.00  0 

-0.996 

C  .0 

65. 

2.522 

30.000 

-1.365 

0.0 

70. 

2.380 

90.000 

-1.867 

CO 

75. 

2.205 

90.000 

-2.531 

o.c 

80. 

1.999 

90.000 

-3.384 

0  . 0 

35. 

1.768 

90.000 

-4.448 

0.0 

90. 

1.524 

90.000 

-5.742 

0.0 

95. 

1.277 

90.000 

-7.275 

r    p 

100. 

1.04  2 

90.000 

-9.045 

0.0 

10  5. 

C.823 

90.000 

-11.030 

0.0 

110. 

0.647 

90.000 

-13.177 

0.0 

115. 

0.5C2 

90.000 

-15.388 

0.0 

120. 

C.393 

90.000 

-17.512 

o.c 

125. 

C.318 

90.000 

-19.351 

o.c 

130. 

C.272 

90.000 

-20.720 

u  • 

13  5. 

C24  8 

9C.000 

-21.525 

0 . 0 

140. 

0.239 

90.000 

-21.817 

o.c 

145. 

C?41 

90.000 

-21.749 

0.0 

150. 

0.249 

90.000 

-21.488 

o.c 

155. 

C.258 

90.000 

-21.164 

L  .  U 

160. 

0.267 

90.000 

-20.853 

0  . 0 

165. 

C.275 

90.000 

-20.597 

0  . 0 

170. 

C.281 

90.000 

-?0.412 

0.0 

175. 

0.285 

90.000 

-20. 3^1 

-V   •   «J 

180. 

0.286 

90.000 

-20.264 

01REOT  IVITY/ISO      FI  3DB      ELEV  30B 
5.77208        78.  35. 

This  DIRECTIVITY  FIGURE  WAS  COMPUTED 

ASSUMING  THAT  THE  PATTERN  IS  WELL-BEHAVED. 
IT  HAS  \C  "FANING  IF  THERE  ARE  SIDE-LOBES. 
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CURRENT    C  IS TR I  PUT ITN     1 

FREQ  = 

0.170E     10 

ELFV 

FI 

E 

PSI 

EOB 

O.G 

CO 

2.330 

90.000 

0.0 

1C. 

CO 

2.257 

9C.000 

-0.275 

20. 

CO 

2.052 

90.000 

-1.103 

30. 

0 . 0 

1.749 

9C.000 

-2.493 

40. 

o.c 

1.394 

9C.0O0 

-4.462 

50. 

CO 

1.035 

90.000 

-7.051 

60. 

o.c 

C708 

90.000 

-10.344 

70. 

CO 

C437 

90.000 

-14.533 

80. 

o.c 

C.233 

9C.000 

-20.013 

90. 

180. 

0.C97 

-90.000 

-27.627 

80. 

ifiO. 

0.025 

-90.000 

-39.335 

70. 

18C. 

0.0C5 

-89.999 

-53.477 

60  . 

180. 

C  .  C  1 6 

-90. 000 

-43.384 

50. 

180. 

CC35 

-90.000 

-36.506 

40. 

18C 

0.045 

-30.000 

-34.311 

30. 

180. 

C.C40 

-90.000 

-35.350 

20. 

18C 

C.024 

-90.000 

-39.586 

10. 

ieo. 

CC0  9 

-9C.000 

-48.263 

0  . 0 

180. 

C.0C3 

-90.000 

-58.853 

o.c 

c.c 

2.330 

9C.000 

0.0 

0.0 

5. 

2.322 

90.000 

-0.029 

0.0 

10. 

2.299 

90.000 

-0.  118 

15. 

2.259 

90.00  0 

-0.259 

0.0 

20, 

2.2C3 

90.000 

-0.486 

o.c 

25. 

2.132 

90.000 

-0.773 

0.0 

30. 

2.043 

9C.000 

-1. 140 

0.0 

35. 

1.939 

90.000 

-1.595 

o.c 

40. 

1.819 

90.000 

-2. 149 

CO 

45. 

1.685 

90.noo 

-2.817 

0.0 

50. 

1.537 

90.000 

-3.614 

1    h 

55. 

1.378 

90.000 

-4.560 

C  •  0 

6  0. 

1.212 

90.00  0 

-5.630 

0.0 

65. 

1.C40 

9C.000 

-7.005 

0.0 

70. 

C.868 

9C.000 

-8.575 

0.0 

75. 

C.7C0 

90.000 

-10.445 

CO 

80. 

C.540 

90.00  0 

-12.699 

CO 

85. 

C.392 

9C.000 

-15.474 

CO 

90. 

C.261 

90.000 

-19.029 

CO 

95. 

C.147 

90.000 

-23.981 

ICO. 

CC54 

90.000 

-32.655 

o.c 

1C5. 

C.C18 

-90.001 

-42. 158 

o.c 

no. 

0.071 

-9C.000 

-30.357 

C  .  0 

115. 

0.1C5 

-90.0QO 

-26.922 

0.0 

120. 

C.123 

-90.000 

-25.516 

125. 

0.129 

-90.000 

-25.138 

CO 

no. 

0.  125 

-90.000 

-25.442 

o.c 

135. 

0.113 

-90.000 

-26.278 

CO 

140. 

0.C97 

-90.000 

-27.579 

o.c 

14  5. 

0.C80 

-90.000 

-29. 327 

o.c 

150. 

CC6  2 

-90.000 

-31.547 

o.c 

155. 

C.045 

-90.000 

-34.306 

G  .0 

16C 

C.G30 

-9C.000 

-37. 736 

CO 

165. 

0.018 

-9C.C00 

-42.058 

o.c 

170. 

C.C10 

-90.000 

-47.606 

O.G 

175. 

C.C04 

-90.000 

-54.423 

o.c 

180. 

C.C03 

-9C.000 

-58.852 

DIRECTIVITY/ISC      FI  3DB      ELFV  3DB 
7.881CP         48.  35. 

THIS  CIPFCTIVITY  FIGURE  WAS  COMPUTED 

ASSUMING  THAT  THE  PATTERN  IS  WELL-BEHAVED. 
IT  HAS  NG  MEAN  INC  IF  THERE  ARE  SIDE-LOBES. 
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1  or 


CURRENT    DISTRIBUTION    2 

FREQ  = 

0.170E    10 

ELEV 

FI 

E 

PSI 

EDB 

O.C 

CO 

2.235 

90.000 

0.0 

10. 

O.C 

2.176 

90.000 

-0.233 

20. 

O.C 

2.CC9 

90.000 

-0.927 

30. 

0.0 

1.760 

90.000 

-2.079 

40. 

CO 

1.461 

90.000 

-3.693 

50. 

CO 

1.144 

90.000 

-5.820 

60. 

CO 

0.831 

90.000 

-8.598 

70. 

CO 

0.541 

9C.000 

-12.327 

8C. 

O.C 

C.295 

9C.OO0 

-17.600 

90. 

leo. 

0.116 

-90.000 

-25.686 

80. 

180. 

C023 

-90.000 

-39.767 

70. 

180. 

0.C14 

-9C.000 

-43.819 

60. 

180. 

C.C66 

-90.000 

-30.559 

50. 

180. 

C.140 

-90.000 

-24.072 

40. 

180. 

C.2CC 

-90.000 

-20.971 

30. 

180. 

0.229 

-90.000 

-19.789 

20. 

180. 

0.231 

-90.000 

-19.724 

1C. 

180. 

C.221 

-90.000 

-20.101 

CO 

180. 

0.216 

-9C.000 

-20.317 

O.C 

CO 

2.235 

9C.000 

0.0 

0.0 

5. 

2.231 

90.000 

-0.019 

0.0 

10. 

2.216 

90.000 

-0.075 

CO 

15. 

2.192 

90.000 

-0. 171 

0.0 

20. 

2.  157 

90.000 

-0.311 

O.C 

25. 

2.111 

90.000 

-0.499 

O.G 

30. 

2.052 

90.000 

-0.742 

O.C 

35. 

1.981 

90.000 

-1.049 

0.0 

40. 

1.895 

90.000 

-1.432 

c.c 

45. 

1.795 

90.000 

-1.907 

0.0 

50. 

1.678 

90.000 

-2.490 

O.C 

55. 

1.545 

9C.C00 

-3.206 

0.0 

60. 

1.397 

90.000 

-4.081 

CO 

65. 

1.235 

9C.000 

-5.153 

0.0 

70. 

1.C61 

90.000 

-6.470 

CO 

75. 

C.880 

90.00  0 

-8.099 

CO 

80. 

C.695 

9C.00C 

-10. 144 

CO 

85. 

C.513 

90.000 

-12.787 

O.C 

90. 

0.338 

90.000 

-16.404 

CO 

95. 

C.177 

90.000 

-22.045 

0    .    v> 

IOC. 

C.C33 

90.000 

-36.655 

0.0 

105. 

C.C90 

-90.000 

-27.934 

0.0 

110. 

C.189 

-90.000 

-21.466 

0.0 

115. 

0.264 

-90.000 

-18.554 

CO 

120. 

0.316 

-90.000 

-16.989 

CO 

125. 

0.347 

-90.000 

-16. 173 

0.0 

130. 

C.360 

-90.000 

-15.850 

CO 

135. 

C.359 

-90.000 

-15.880 

0.0 

140. 

C.347 

-90.000 

-16.172 

CO 

145. 

C.328 

-90.000 

-16.657 

CO 

150. 

C.306 

-90.000 

-17.276 

0.0 

155. 

C.282 

-90.000 

-17.968 

CO 

160. 

C.260 

-90.000 

-18.672 

0.0 

165. 

0.242 

-90.000 

-19.322 

c.c 

170. 

0.227 

-90.000 

-19.850 

CO 

175. 

C.219 

-90.000 

-20. 196 

CO 

180. 

0.216 

-90.000 

-20.317 

DIRECTIVITY/ISO      FI  3DB      ELEV  30B 
7.4510B        53.         35. 

THIS  DIRECTIVITY  FIGURE  WAS  COMPUTED 

ASSUMING  THAT  THE  PATTERN  IS  WELL-BEHAVED. 
IT  HAS  NO  MEANING  IF  THERE  ARE  SIDE-LOBES. 
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Toe 


CURRENT  DISTRIBUTE  1 

FRE0= 

0.200E  10 

ELFV 

FI 

E 

PSI 

ED 

O.C 

CO 

2, 

,581 

-9C. 

,000 

-0. 

000 

10. 

O.C 

2. 

,494 

-90, 

,000 

-0. 

297 

20. 

o.c 

2. 

,249 

-90. 

,000 

-1. 

196 

30. 

CO 

1, 

,889 

-90. 

,000 

-2. 

710 

40. 

CO 

1. 

,475 

-90. 

,00  0 

-4, 

861 

50. 

CO 

1. 

,067 

-9C. 

,  0n0 

-7. 

672 

6C. 

c.c 

C 

,714 

-90, 

,000 

-11. 

161 

70. 

CO 

0, 

,440 

-90, 

,00  0 

-15. 

363 

80. 

CO 

o. 

,247 

-90, 

,000 

-20. 

384 

90. 

180. 

c. 

,  120 

90. 

,000 

-26. 

630 

80. 

180. 

C 

,040 

90. 

,000 

-36. 

114 

7C. 

180. 

C 

,C12 

-90. 

,001 

-46. 

740 

60. 

180. 

C. 

,052 

-90. 

,00  0 

-33. 

987 

50. 

180. 

0, 

,C89 

-90. 

,000 

-29. 

259 

40. 

18C 

C 

,128 

-9C 

,000 

-26. 

067 

30. 

180. 

0. 

,169 

-90. 

,000 

-23. 

684 

20. 

180. 

C. 

,2C5 

-9C. 

,000 

-21. 

993 

1C. 

180. 

c. 

,231 

-9  0, 

,000 

-20. 

982 

0.0 

18C. 

c, 

,24C 

-9C. 

,000 

-20. 

645 

0.0 

CO 

2. 

,581 

-90, 

,Ono 

0. 

0 

0.0 

5. 

2, 

,572 

-90. 

,00  0 

-0. 

032 

0.0 

10. 

2, 

,543 

-90. 

,000 

-0. 

129 

CO 

15. 

2, 

,496 

-90. 

,00  0 

-0. 

293 

0.0 

20. 

2, 

,429 

-9C. 

,000 

-0. 

527 

P  0 

25. 

2. 

,344 

-90. 

,0^0 

-0. 

836 

CO 

30. 

2, 

,242 

-9C. 

,  OJO 

-1. 

225 

o.c 

35. 

2. 

,  122 

-90. 

,000 

-1. 

702 

o.c 

40. 

1. 

,987 

-9C. 

,000 

-2. 

275 

0.0 

45. 

1. 

,838 

-90. 

,000 

-2. 

952 

0.0 

50. 

1. 

,677 

-9C. 

,000 

-3. 

743 

CO 

55. 

1. 

,509 

-90. 

,00  0 

-4. 

660 

CO 

60. 

1, 

,337 

-90. 

,00  0 

-5. 

712 

CO 

65. 

1. 

,165 

-90. 

,000 

-6. 

910 

CO 

70. 

C 

,997 

-90, 

,000 

-8. 

263 

CO 

75. 

C 

,838 

-9C, 

,000 

-9. 

776 

0.0 

CO  » 

0. 

,691 

-90, 

,000 

-11. 

450 

CO 

85. 

c 

,560 

-90. 

,000 

-13. 

276 

O.C 

90. 

c< 

,447 

-90, 

,000 

-15. 

226 

CO 

95. 

c 

,354 

-90, 

,000 

-17. 

248 

0.0 

ICC 

c. 

,281 

-9C, 

,00  0 

-19. 

249 

o.c 

105. 

c. 

,228 

-90, 

,000 

-21. 

092 

CO 

110. 

c, 

,191 

-90. 

,coo 

-22. 

608 

CO 

115. 

c 

,170 

-90, 

,000 

-23. 

649 

CO 

120. 

c. 

,  160 

-9C, 

,000 

-24. 

153 

CO 

125. 

c 

,160 

-90, 

,000 

-24, 

178 

CO 

130. 

c, 

,165 

-90, 

,000 

-23. 

867 

CO 

135. 

0, 

,175 

-90. 

,000 

-23. 

375 

o.c 

140. 

0. 

,  186 

-90, 

,000 

-22. 

826 

CO 

145. 

c. 

,198 

-90, 

,000 

-22. 

298 

o.c 

150. 

0. 

2C9 

-9C, 

,coo 

-21. 

832 

0.0 

155. 

c 

,219 

-90, 

,000 

-21. 

446 

o.c 

160. 

0. 

,226 

-9C, 

.000 

-21. 

143 

CO 

165. 

c 

,232 

-90. 

,000 

-20. 

917 

CO 

170. 

c. 

236 

-90. 

,000 

-20. 

764 

CO 

175. 

c 

239 

-90. 

,000 

-20. 

674 

0.0 

180. 

0. 

240 

-90. 

,000 

-20. 

645 

DIRECTIVITY/ISO      FI  3CB      ELEV  3DB 
7.881DP        48.         35. 

THIS  DIRECTIVITY  FIGURE  MS  COMPUTED 

ASSUMING  THAT  THE  PATTERN  IS  WELL-BEHAVED, 
IT  HAS  NO  MEANING  IF  THERE  ARF  SIDE-LOBES, 


98 


CURRENT    DISTKlBliT  K'N    2 

FREQ  = 

0.2C0E    1^ 

tLEV 

FI 

E 

PSI 

EDB 

G  .  C 

0.0 

2.2P8 

-Qp.000 

-0.000 

10. 

CO 

2.201 

-90.000 

-0. 335 

2C 

o.c 

1.959 

-90.000 

-1.351 

30. 

c.c 

1.607 

-90.000 

-3.071 

40. 

0.0 

1.210 

-90.000 

-5.530 

50. 

CO 

0.835 

-90.000 

-8. 757 

60. 

l>      r 

0.527 

-9C.000 

-12.745 

70. 

o.c 

0.308 

-90.000 

-17.409 

P  ). 

CO 

0.171 

-90.000 

-22.546 

SO. 

IPC 

CC90 

90.000 

-28.079 

80. 

IP  p. 

CC3P 

90.000 

-35.573 

70. 

180. 

C.010 

-90.001 

-47. 134 

60. 

180. 

CC69 

-90.000 

-30.415 

50. 

IPC 

0.143 

-90.000 

-24.070 

40. 

IPC 

C22P 

-9C.0O0 

-20.013 

30. 

IPC 

0.314 

-90.000 

-17.245 

20. 

IRC 

0.3PB 

-90.000 

-15.423 

10. 

IPO. 

0  .437 

-90.00  0 

-14. 385 

C.P 

IPC 

C454 

-90.00  0 

-14.048 

C  .  0 

0.0 

2.2P8 

-90.000 

0.0 

0  •  c 

5. 

2.273 

-90.000 

-0.040 

0.0 

10. 

2.246 

-90.000 

-0.160 

0 .  c 

15. 

2.  194 

-90.000 

-0.363 

0.0 

20. 

2.123 

-90.000 

-0.651 

0.0 

25. 

2.032 

-90.000 

-1.031 

0.0 

30. 

1.923 

-90.000 

-1.509 

CO 

35. 

1.798 

-90.000 

-2.093 

0.0 

AC 

1.659 

-90.000 

-2. 791 

0.0 

45. 

1.5G9 

-90.000 

-3.615 

0.0 

50. 

1.351 

-90.000 

-4.574 

c.c 

55. 

1.190 

-90.000 

-5.682 

0.0 

60. 

1.C28 

-90.000 

-6.947 

V  •  o 

65. 

0.872 

-90.000 

-8.380 

0.0 

70. 

0.725 

-9C.000 

-9.984 

CO 

75. 

0.592 

-90.000 

-11. 750 

c.c 

PC 

C.475 

-90.000 

-13.651 

0  .0 

85. 

0.379 

-90.000 

-15.625 

CO 

90. 

0.30  3 

-90.000 

-17.555 

c.c 

95. 

C249 

-90.000 

-19.257 

o.c 

100. 

0.216 

-90.000 

-20.506 

^  •  c 

105. 

C.2C1 

-90.000 

-21.124 

c.c 

11C 

C.2C2 

-90.000 

-21.084 

CO 

1  15. 

C.215 

-9C.000 

-20.528 

C  .0 

120. 

0.238 

-90.000 

-19.674 

CO 

125. 

0.26  5 

-90.000 

-18. 715 

c.c 

130. 

0.296 

-90.000 

-17.778 

CO 

135. 

C.326 

-90.000 

-16.933 

c.c 

14C. 

0.354 

-90.000 

-16.208 

c.c 

145. 

0.379 

-90.000 

-15.610 

o.c 

15C 

0.401 

-90.000 

-15. 132 

CO 

155. 

C.418 

-90.000 

-14.760 

CO 

16C 

0.432 

-90.000 

-14.482 

c.c 

165. 

0.442 

-90.000 

-14.282 

CO 

170. 

C.449 

-90.000 

-14. 148 

o.c 

175. 

C453 

-90.000 

-14.072 

0  .0 

IPC 

C.454 

-90.000 

-14.048 

JIRECTIVITY/ISO      Fl  ^DB      ELEV  3DB 
9.820CB        43.  25. 

THIS  CIRECTIVITY  FIGURE  WAS  COMPUTED 

ASSUMING  THAT  THt  PATTERN  IS  WELL-BEHAVED. 
IT  HAS  NO  MEANING  IF  THERE  ARE  SIDE-LOBES. 
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CURRENT  DISTRIBUTION  I 

FREQ  = 

0.230E  10 

ELEV 

FI 

E 

PSI 

ED 

0.0 

0.0 

2, 

,470 

90. 

.000 

-0. 

000 

10. 

0.0 

2, 

,400 

90. 

,00  0 

-0. 

250 

20. 

0.0 

2, 

,20  1 

90, 

.000 

-1. 

002 

30. 

CO 

1. 

,9C3 

9C. 

.000 

-2. 

266 

40. 

0.0 

1. 

,547 

90. 

,000 

-4. 

064 

50. 

CO 

1, 

,178 

90, 

,000 

-6. 

435 

60. 

CO 

C, 

,831 

90, 

.000 

-9. 

467 

70. 

CO 

c 

,532 

90, 

.000 

-13. 

341 

ac. 

o.O 

0, 

,295 

9C. 

,000 

-18. 

447 

90. 

180. 

0. 

,128 

-90. 

,000 

-25. 

722 

80. 

180. 

c. 

,C29 

-89, 

,999 

-38. 

723 

70. 

IPC 

c. 

,010 

89. 

,999 

-48. 

285 

60. 

180. 

o, 

,001 

90, 

,000 

-72. 

899 

50. 

180. 

c. 

,037 

-90. 

,000 

-3  6. 

419 

40. 

180. 

c. 

,086 

-90, 

,000 

-29. 

140 

30. 

180. 

c. 

,133 

-9C, 

.000 

-25. 

396 

20. 

180. 

0. 

,  169 

-9C. 

.000 

-23. 

306 

10. 

18C. 

0, 

,191 

-90. 

,000 

-22. 

228 

o.c 

180. 

c 

,  199 

-90, 

,00  0 

-21. 

894 

CO 

O.C 

2. 

,470 

90. 

,00  0 

0. 

0 

CO 

8. 

2. 

463 

90. 

,coo 

-0. 

027 

CO 

10. 

2. 

,440 

90. 

,000 

-0. 

107 

0.0 

18. 

2. 

,402 

90. 

,000 

-0. 

244 

O.C 

20. 

2. 

,348 

90. 

,000 

-0. 

441 

CO 

25. 

2. 

,278 

90. 

,000 

-0. 

704 

CO 

50. 

2. 

,  192 

90, 

,00^ 

-1. 

039 

CO 

35. 

2. 

,089 

90. 

.000 

-1. 

,454 

0  .  G 

4C 

1. 

,971 

90. 

,00  0 

-1. 

960 

CO 

45. 

1, 

,838 

90. 

,000 

-2. 

567 

CO 

50. 

1, 

,692 

90. 

,000 

-3. 

288 

0.0 

55. 

1. 

,535 

90, 

,000 

-4. 

135 

CO 

60. 

1, 

,369 

90, 

,000 

-5, 

124 

w  .  U 

65. 

1, 

,200 

90. 

,000 

-6. 

274 

Q.O 

70. 

1. 

,029 

90, 

,000 

-7. 

606 

c.c 

75. 

0. 

,862 

90, 

,000 

-9. 

,  145 

O.C 

80. 

C. 

702 

90, 

,000 

-10. 

928 

CO 

85. 

0, 

,553 

90. 

,000 

-13. 

005 

0.0 

90. 

c 

417 

90. 

,000 

-15. 

455 

CO 

95. 

c. 

296 

9C, 

,000 

-18. 

418 

CO 

100. 

0. 

192 

90, 

,000 

-22. 

178 

CO 

105. 

c 

,  105 

90, 

,000 

-27. 

470 

0.0 

110. 

c 

,032 

90. 

,000 

-37. 

619 

Q.O 

115. 

c 

025 

-90. 

,000 

-39. 

817 

O.C 

120. 

C. 

C70 

-90, 

,00  0 

-30. 

901 

Q.C 

125. 

c 

1C5 

-90. 

,000 

-27. 

426 

CO 

130. 

C. 

131 

-90. 

.000 

-25. 

500 

CO 

135. 

c 

150 

-90. 

,000 

-24. 

305 

CO 

14C 

0. 

,165 

-90. 

,000 

-23. 

520 

O.C 

145. 

C. 

175 

-90, 

,000 

-22. 

988 

0.0 

150. 

G. 

183 

-90. 

,000 

-22. 

618 

O.C 

155. 

c 

188 

-90, 

,000 

-22. 

358 

c.c 

160. 

c 

192 

-90. 

,000 

-22. 

173 

0.0 

165. 

c. 

195 

-90. 

,000 

-22. 

044 

0.0 

170. 

0. 

197 

-90, 

,000 

-21. 

959 

c.c 

175. 

c. 

198 

-90, 

,000 

-21. 

910 

O.C 

180. 

c. 

199 

-90, 

,000 

-21. 

894 

OIRECTIVITY/tSO      FI  3CB      ELEV  3DB 
7.881De        48.         35. 

THIS  DIRECTIVITY  FIGURE  WAS  COMPUTED 

ASSUMING  THAT  THE  PATTERN  IS  WELL-BEHAVED. 
IT  HAS  NG  MEANING  IF  THERE  ARE  SIDE-LOBES. 
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fo? 


CURRENT    CISTRIRUTICN    2 

FREQ  = 

0.230E    10 

ELEV 

FI 

E 

PSI 

EDB 

0.0 

CO 

2.470 

90.000 

-0.000 

1C. 

CO 

2.400 

90.000 

-0.250 

20. 

c.c 

2.201 

90.000 

-1.002 

30. 

c.c 

1.903 

9C.000 

-2.266 

40. 

c.c 

1.547 

90.000 

-4.064 

50. 

CO 

1.178 

90.000 

-6.435 

60. 

0.0 

0.831 

9C.000 

-9.467 

7C. 

CO 

0.532 

90.000 

-13.341 

80. 

CO 

0.295 

90.000 

-18.447 

90. 

180. 

C128 

-90.00C 

-25.722 

80. 

180. 

CC29 

-89.999 

-38.723 

70. 

180. 

0.010 

89.999 

-48.285 

60. 

180. 

0.CC1 

90.000 

-72.899 

50. 

180. 

C037 

-90.000 

-36.419 

40. 

180. 

C.C86 

-90.000 

-29. 140 

30. 

190. 

0.  133 

-90.000 

-25.396 

20. 

180. 

C.169 

-90.000 

-23.306 

10. 

180. 

0.191 

-90.000 

-22.228 

0.0 

lao. 

0.199 

-90.000 

-21.894 

C.C 

c.c 

2.470 

9C.O00 

0.0 

0.0 

5. 

2.463 

9C.000 

-0.027 

c.c 

10. 

2.440 

90.000 

-0. 107 

CO 

15. 

2.402 

90.000 

-0.244 

0.0 

20. 

2.348 

90.000 

-0.441 

0.0 

25. 

2.278 

90.000 

-0.704 

0.0 

30. 

2.  192 

90.000 

-1.0  39 

CO 

35. 

2.C89 

90.000 

-1.454 

0.0 

40. 

1.971 

9C.00C 

-1.960 

0.0 

45. 

1.838 

9C.C0O 

-2.567 

CO 

50. 

1.692 

90.000 

-3.288 

CO 

55. 

1.535 

9  0.000 

-4.  135 

cc 

60. 

1.369 

90.000 

-5. 124 

0.0 

65. 

1.200 

9C.00C 

-6.274 

0.0 

70. 

1.C29 

90.000 

-7.606 

0.0 

75. 

0.862 

90.000 

-9.145 

o.c 

80. 

0.702 

90.000 

-10.928 

0 . 0 

85. 

C.553 

9C.00  0 

-13.005 

CO 

90. 

0.417 

90.000 

-15.455 

CO 

95. 

0.296 

9C.000 

-18.418 

o.c 

100. 

0.192 

90.000 

-22.178 

0     ^ 

105. 

C.105 

90.000 

-27.470 

CO 

110. 

0.032 

90.000 

-37.619 

0.0 

115. 

0.025 

-90.000 

-39.817 

0.0 

120. 

C.07C 

-90.010 

-30.901 

CO 

125. 

C.105 

-90.000 

-27.426 

0.0 

130. 

0.131 

-90.000 

-25.500 

0.0 

135. 

C.  150 

-90.000 

-24.305 

CO 

140. 

0.165 

-90.000 

-23.520 

0.0 

145. 

C175 

-90.000 

-22.988 

o.c 

150. 

C.183 

-90.000 

-22.618 

o.c 

155. 

0.188 

-90.000 

-22.358 

o.c 

16C 

C192 

-90.000 

-22.173 

CO 

165. 

0.  195 

-90.000 

-22.044 

o.c 

170. 

C197 

-90.000 

-21.959 

0.0 

175. 

C198 

-90.000 

-21.910 

o.c 

180. 

0.199 

-90.000 

-21.894 

CIRECTIVITY/ISO      FI  3CB      ELEV  30B 
7.881DB        48.  35. 

THIS  DIRECTIVITY  FIGURE  fcAS  COMPUTED 

ASSUMING  THAT  THE  PATTERN  IS  kELL-REHA VEQ. 
IT  HAS  NC  MEANING  IF  THERE  ARE  SIDE-LORES. 
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II.   Program  SMIGUEL 

This  program  also  is  written  in  code  FORTRAN  IV  for  the  IBM  360 
system. 

The  program  uses  seven  subroutines,  as  explained  in  SAMPLE  PROGRAM 
SMIGUEL  that  follows.   The  subroutines  supply  to  the  main  program  the 
characteristics  of  the  line  and  the  non-uniformity  corrections.   The  main 
program  uses  equations  33,  first  to  find  the  input  impedance,  and  then  to 
find  the  current  distribution. 

For  each  set  of  input  data  the  impedance  and  current  distribution  are 
computed  twice — one  time  without  accounting  for  mutual  impedance  effects, 
and  a  second  time  accounting  for  them. 

Again,  the  program  has  a  large  number  of  comment  cards  that,  hope- 
fully, make  it  self-explanatory. 
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